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Cold Rydberg atoms exposed to strong magnetic fields possess unique properties which open the 
pathway for an intriguing many-body dynamics taking place in Rydberg gases consisting of either 
matter or anti-matter systems. We review both the foundations and recent developments of the 
field in the cold and ultracold regime where trapping and cooling of Rydberg atoms have become 
possible. Exotic states of moving Rydberg atoms such as giant dipole states are discussed in detail, 
including their formation mechanisms in a strongly magnetized cold plasma. Inhomogeneous field 
configurations influence the electronic structure of Rydberg atoms, and we describe the utility of 
corresponding effects for achieving tightly trapped ultracold Rydberg atoms. We review recent 
work on large, extended cold Rydberg gases in magnetic fields and their formation in strongly 
magnetized ultracold plasmas through collisional recombination. Implications of these results for 
current antihydrogen production experiments are pointed out, and techniques for trapping and 
cooling of such atoms are investigated. 
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I. INTRODUCTION 

Rydberg atoms, by virtue of their special properties, 
are of fascination in atomic and optical physics [l[ . These 
atoms whose size scales as the square of the Rydberg 
principal quantum number n, are ideally suited for fun- 
damental quantum interrogations as well as detailed clas- 
sical analyses which test the correspondence principle. 

Owing to the large excursions of the bound Rydberg 
electron, highly excited states have extremely long life- 
times, scaling as n 3 . Their exaggerated dimensions also 
imply large geometrical cross sections, such that atomic 
processes like Penning ionization or inelastic electron 
collisions (ionization, recombination and (de)-excitation) 
occur with high rates in an ultracold plasma environ- 
ment. Moreover, the large polarizabilities of Rydberg 
atoms makes them highly susceptible to external elec- 
tric and magnetic fields — the latter being the primary 
subject of the present article. 

The largest discovered Rydberg atoms present in a nat- 
ural environment were carbon atoms in interstellar space 
with a principal quantum number of n ~ 1009 2] which 
have been detected due to their radio recombination lines. 
Rydberg atoms in cool tenuous interstellar media offer a 
valuable opportunity for measuring magnetic fields. Mea- 
surements of Hna lines (An = 1 radio recombination 
transitions) allow for the determination of the Zeeman 
splitting in hydrogen, from which typically small mag- 
netic field strengths in the range of fiG can be inferred 
HU. Although small in strength, such magnetic fields 
in tenuous environments, where neutral densities are less 
than 10 4 cm -3 , can play an important role for the inter- 
stellar gas dynamics Q. 

The laboratory creation of cold and dilute highly- 
excited media has become possible with the develop- 
ment of techniques to laser-cool and trap neutral atoms 
@; S 0, H ■ Laser-excitation of cold atomic gases leads 
to an intriguing many-body system with very strong long- 
range interactions. Together with their extremely long 
lifetimes the strong interactions between Rydberg atoms 
make for promising candidates in quantum information 
processing [l(| EH EH EH EH- A range of questions, 
arising from the collective excitation dynamics in cold 
Rydberg gases 0, EI El EE H M M M M H3 and 
the formation and non-e quili brium evolution of ultracold 
neutral plasmas pj, M,M, El, El, M, HI, EH, H, E3 
have been studied in experiment and theory. Since out- 
lining these emerging topics is certainly beyond the scope 
of the present article, we refer to some recent review ar- 
ticles [25j, |26l [27J for more comprehensive discussions. 

A very different class of Rydberg atoms yet appears in 
strong magnetic fields, as realized through recombination 
processes in cold antimatter plasmas at CERN [28|, [29| , 
in ion storage rings [3(j, and in cold atom traps, sub- 
ject to magnetic fields of several Tesla [3l|. Due to their 
fragility, the effects of such laboratory magnetic fields 
on the electronic structure of Rydberg atoms are simi- 
larly pronounced as those of very strong, extraterrestrial 



fields (~ 10 4 Tesla) [32| on ground state atoms. Highly 
excited Rydberg states, hence, provide a well-accessible 
venue for combined theoretical and experimental studies 
of strong magnetic field effects on atomic systems. From 
a different perspective, highly inhomogeneous magnetic 
fields may offer a promising possibility for a controlled 
manipulation of Rydberg atoms, which are otherwise dif- 
ficult to trap by standard optical techniques, developed 
for ground state atoms. 

Here, we will focus on magnetic field induced static 
and dynamical effects on Rydberg atoms, i.e. on their 
internal and translation degrees of freedom, their forma- 
tion in cold plasmas as well as their evolution in external 
traps. This review article is organized as follows. Section 
Hfl is devoted to existing experiments in magnetized Ry- 
dberg gases and plasmas, with emphasis on the atomic 
excitation from laser-cooled gases of ground state atoms. 
Antihydrogen production in Penning traps through re- 
combination in magnetized plasmas is briefly outlined 
here. Section lllll first focuses on the fundamental prop- 
erties of atomic systems in homogeneous external mag- 
netic fields. Following a discussion of the symmetries and 
constants of motion, the latter are exploited to perform 
a gauge-independent pseudoseparation of the center of 
mass motion of neutral atoms in magnetic fields or com- 
bined perpendicular electric and magnetic fields. Effects 
due to the magnetic field-induced coupling of the center 
of mass motion to the internal motion are reviewed on a 
classical-mechanics level. The main focus of this section 
are the giant dipole states whose preparation, properties 
as well as applications to matter-antimatter systems are 
studied in detail. We conclude with a discussion of the 
concept of the guiding center approximation to describe 
strongly magnetized atoms. In section IIVI the structure 
and quantum dynamics of atomic systems in inhomo- 
geneous magnetic fields is addressed. We first describe 
their properties assuming immobile nuclei, followed by 
a systematic approach to trapping of ultracold Rydberg 
atoms in 3D quadrupole and Ioffe-Pritchard traps. An 
analysis of the lifetime of trapped Rydberg atoms and 
their long-time dynamics in strong field traps comple- 
ment this section. Section [V] is dedicated to Rydberg 
atom formation and subsequent collisions in cold magne- 
tized plasmas, with particular reference to antihydrogen 
experiments. We conclude with an outlook in section [VTl 
pointing out some perspectives for future studies. 



II. MAGNETIZED RYDBERG GASES AND 
PLASMAS 



Over the last decade, experimental efforts in pho- 
toexcitation and photoionization of laser-cooled atoms 
have led to creation of cold highly-excited Rydberg gases 
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34 and ultracold neutral plasmas 
4ll |42| . At strong magnetic fields 
B and sufficiently low thermal velocities Vth — \/k^T '/m 
a situation can be realized, for which the cyclotron ra- 
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dius r c = uth/^c becomes much smaller than the dis- 
tance of closest approach b — e 2 /k B T, and the cyclotron 
frequency lo c = eB/mc provides the shortest dynamical 
times scale of charges with mass m. Such strongly mag- 
netized conditions are naturally realized in non-neutral 
plasmas of either electrons or laser-cooled ions, confined 
in the strong magnetic fields of a Penning trap [43| . 
Their intense study over the last few decades has led, 
for example, to a deeper understanding of strongly cor- 
related plasmas and various transport processes. Neu- 
trality, as realized in cold Rydberg gases, adds another 
degree of complexity, by providing a collisionally rich en- 
vironment in which phenomena such as an avalanche- like 
plasma formation [44J, angular momentum mixing col- 
lisions [45[, Penning ionization and Rydberg atom 
formation through three-body recombination [421 ] have 
been observed. 



ries of sharp resonances, at frequency intervals of 3w c /2 
and extending well into the zero-field ionization contin- 
uum. In contrast regular Landau levels are spaced by the 
electron cyclotron frequency lu c = eB/mc. The contin- 
uum features in the spectra of magnetized atoms — later 
termed quasi-Landau resonances — have been studied in 
a number of subsequent experiments (5ll . l52l . |53| . The 
diamagnetic Rydberg atom in strong magnetic fields has, 
moreover, received considerable attention as a paradigm 
system for quantum chaos (54|, permitting detailed and 
precise comparisons between theory and experiment [55| . 
For further details on these aspects of Rydberg atoms in 
strong magnetic fields we refer the reader to the corre- 
sponding reviews [S(| H3, HH and cited references 
therein. 



B. Cold gases in strong magnetic fields 



A. Atomic excitation in strong magnetic fields 

Rydberg atoms are generally produced by collision 
(charge particle excitation and electron exchange) or op- 
tical excitation. The production of cold Rydberg atoms is 
typically achieved through optical excitation. Here laser 
cooled gases in a magneto-optical [H, [TH, [I?], El] or a 
dipole [l8| | trap are used as a starting point, which sub- 
sequently are excited to Rydberg states either by UV ex- 
citation or by two-photon excitation [l(| [fT], [H, Oil • 
Typically two-photon excitation of Rubidium atoms pro- 
ceeds via excitation of the 5iSi/2 ground state to the first 
excited 5P3/2, which is coupled by a second laser to ei- 
ther nSi/2 or nD 3 /2 : 5/2 Rydberg states. The maximum 
excitation efficiency achieved with such laser-pulses is, 
however, typically limited by a high sensitivity to spa- 
tial laser-pulse inhomogeneities and fluctuations in laser 
power and frequency. Potentially much higher efficiency 
can be achieved using the method of stimulated Raman 
adiabatic passage (STIRAP), in which the upper transi- 
tion pulse is applied ahead of the lower one, in an coun- 
terintuitive pulse sequence [47| . The method utilizes the 
existence of a dark state in a three-level system, which is 
composed only of the ground and excited Rydberg state. 
By adiabatically following this dark state during the evo- 
lution of the laser pulse, in principle 100% excitation 
efficiency can be achieved under perfect adiabatic con- 
ditions. In cold Rydberg gases this approach was first 
demonstrated by Cubel et al. [48} . who reported an ex- 
citation efficiency of ~ 70%. 

The effects of strong magnetic fields, i.e. of the dia- 
magnetic coupling, on the electronic structure of Ryd berg 
atoms were originally studied by Jenkins and Segre [491 ] 
through absorption spectroscopy of sodium and potas- 
sium atoms for principal quantum numbers of n < 35. 
A later refinement and extension to Rydberg states of 
up to n = 75 by Garton and Tomkins [50j revealed im- 
portant additional features in the absorption spectrum 
of the Barium I series. Unexpectedly, they found a se- 



Magnetic trapping of ground state atoms, as an impor- 
tant prerequisite to forming magnetically trapped Ryd- 
berg gases, has been demonstrated in the group of Raithel 
[26l [59} . Several millions of rubidium atoms were ac- 
cumulated and laser cooled in a superconducting Ioffc- 
Pritchard trap operating at large bias fields of up to 2.9 
Tesla. 

Strongly inhomogeneous magnetic fields whose magni- 
tude varies over length scales comparable to the size of 
a Rydberg atom drastically alter its electronic structure 
and lead to intricate coupling between the atom's inter- 
nal and center-of-mass dynamics. As will be worked out 
in detail in section HVl such configurations provide very 
tight Rydberg atom confinement with laboratory mag- 
netic fields. 

Ref. [3ll | reported Rydberg atom confinement in mag- 
netic traps with weaker field gradients, based on the 
trapping configuration described above [59| . The experi- 
ment could monitor collective atom cloud oscillations in 
the trap, which was used to extract average magnetic 
moments of the trapped Rydberg atoms. Importantly, 
the experiment found extremely long Rydberg gas life- 
times of up to 200 ms. This observation can be traced 
back to a drastically suppressed spontaneous decay rate 
due to the external magnetic field, which will be dis- 
cussed in detail in section IIVBI The combined mea- 
surement of the lifetimes and magnetic moments pro- 
vided strong experimental evidence for the production 
of so-called guiding center atoms. Such atoms are char- 
acterized by large regular electron orbits transverse to 
the magnetic field lines, superimposed by tightly con- 
fined cyclotron oscillations . The basic properties of 
nearly circular — and, hence, long-lived — guiding cen- 
ter atoms and their response to additional electric will 
be worked out in section IIIIFl While the experimen- 
tal two-photon excitation scheme, described in the pre- 
vious section, exclusively produces Rydberg states with 
low angular momenta, the observed high-angular momen- 
tum, circular guiding center are formed through /-mixing 
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collisions with free plasma electrons. There are several 
processes that lead to the formation of free electrons in 
a cold Rydberg gas, whose relative importance is cur- 
rently studied intensively by several grou ps in zero-field 
Rydberg gases M, M, M, H H 1| 

An even more collisionally rich and strongly magne- 
tized environment has been realized in [66| by trapping 
an ultracold neutral rubidium plasma in the strong mag- 
netic fields of a, so-called nested Penning trap (see next 
section). Beyond providing a well suited platform to 
study interesting collective and collisional dynamics in 
the strongly magnetized regime, such experiments can 
also permit important case studies for ongoing antihy- 
drogen experiments, that produce antimatter hydrogen 
atoms through collisional recombination in ultracold two- 
component positron-antiptoton plasmas p8l . |29| . 

C. Antihydrogen production from cold, magnetized 
plasmas 

Several collaborations (ATRAP 0, ATHENA [H, 
ALPHA [13], ASACUSA [H], AEGIS [11]), operating 
at the Antiproton Decelerator at CERN, are working on 
the production of cold antihydrogen atoms, to provide an 
important testing ground for fundamental physics. One 
produced, cold antihydrogen atoms could be employed 
for comparative high precession spectroscopy with the 
hydrogen atom for accurate tests of CPT (charge, parity 
and time) invariance [lO, [7l|, [72|, [73| and could be used to 
measure gravitational forces between neutral matter and 
antimatter, i.e. to test the weak equivalence principle 

The basic ingredient for producing cold antihydro- 
gen (H ) are cold antiprotons (p ) and positrons (e + ). 
The techniques for cooling the initially highly energetic 
antiprotons, produced at CERN, have been developed 
over the two decades by Gabrielse and co-workers (see 
e.g. 73]). In brief, initially, fast p ions are matter 
degraded and captured in a Penning trap and sympa- 
thetically cooled by simultaneously confined electrons, 
which dissipate energy through synchrotron radiation 




FIG. 1: A nested Penning trap. The p accumulate in the 
potential valleys, while the e + are deposited at 4K at the 
saddle of the Penning potential. The magnetic field axis is 
shown. 



[75j |. Positrons are typically obtained from a sodium 
source [79 and later cooled in subsequent steps (see e.g. 
[77L l78l TtOI] for more details). 

In order to simultaneously confine the cold positron 
and antiproton plasma in the same spatial region 
Gabrielse et al. developed a specially designed magnetic 
trap, termed a nested Penning trap [801 ] - The basic prin- 
ciple is schematically shown in FigfTJ As in a standard 
Penning trap radial confinement is provided by a ho- 
mogenous magnetic field oriented along the z-direction 
[43T ]. Axial confinement is achieved through a sequence 
of ring electrodes that realize an effective double well 
potential. Since antiprotons and positrons have oppo- 
site charge, antiprotons are accumulated in the side- wells 
while the positrons are trapped in the central maximum 
of the double well potential. 

Initiating the H formation requires a spatial overlap 
between the two plasma components, which has been re- 
alized experimentally using several different approaches 
[H [H, [MEl. In the first ATHENA and ATRAP ex- 
periments [2et [29j , antiprotons were launched into the 
trapping region with energies above the central well such 
that they would oscillate throughout the trap [H, [29| (see 
FigfTJ) . During each passage through the positron cloud 
antiprotons are cooled via collisions and upon reaching 
sufficiently low velocities start to form highly excited an- 
tihydrogen atoms, predominantly via three-body recom- 
bination (see section [V| . In this approach the velocities 
of formed atoms critically depends on the characteristic 
time of recombination compared to the coo ling timescale, 
and was later found to be rather high [H, l84| . 

In order to achieve lower atomic velocities, a driven 
production scheme was developed in [81j. It starts from 
antiprotons, which are initially cooled into the side-wells 
of the trap, and subsequently heated by rf-fields such that 
they slowly come into overlap with the central positron 
cloud (see Fig.[Tp). In section IVCl we give a further dis- 
cussion of this method and a theoretical interpretation 
of the observed H velocities [85|, |86[ . Finally, in a recent 
experiment the ATRAP collaboration used a third ap- 
proach, based on a gradual increase of the central well 
voltage, which slowly merges both plasma components 
[s3 j. The various detection techniques used to probe the 
properties of p roduced H atoms are reviewed for example 
in [73|,|87|,[If 

To be applicable for the aforementioned spectroscopy 
and gravitation measurements the produced, H atoms 
necessarily have to be translationally cold and in their 
ground state. Since atom formation in cold plasmas pre- 
dominantly proceeds via three-body capture into highly 
excited Rydberg states, the production efficiency of trap- 
pable ground state antihydrogen and their resulting ve- 
locities is crucially determined by the properties of highly 
excited Rydberg atoms and their recombination dynam- 
ics in strong magnetic fields (see section IIVB1 IIV CI and 
|V|) . For a comprehensive overview over the multitude of 
atomic processes occurring in antihydrogen experiments 
we also refer to a recent review by Robicheaux (89j . 
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III. ATOMIC PROPERTIES IN 
HOMOGENEOUS MAGNETIC FIELDS 

The present chapter aims at reviewing the fundamen- 
tal properties of moving atomic systems in homogeneous 
magnetic and/or crossed homogeneous electric and mag- 
netic fields. They are of immediate relevance to the 
dynamics and processes taking place in matter as well 
as antimatter Rydberg gases and plasmas in a strong 
magnetic field and determine many of their properties. 
Firstly, i.e. in section Ull A\ we will introduce the Hamil- 
tonian for both neutral and charged atomic systems and 
discuss its symmetries and resulting constants of motion. 
Interpretations of the conserved pseudomomentum are 
provided. Section IIIIBI is devoted to the pseudosepara- 
tion of the center of mass motion with a focus on neutral 
systems and an outline on ions. The gauge- independent 
pseudoseparation provides us with a generalized poten- 
tial whose properties are discussed firstly here. Dynami- 
cal effects due to the coupling of the center of mass and 
internal (electronic) motion in a homogeneous magnetic 
field are studied in section IIII CI for the underlying clas- 
sical dynamics. The characteristics for both the case of 
a vanishing and a nonvanishing pseudomomentum are 
explored. The giant dipole states emerging due to the 
double- well structure of the generalized potential for suf- 
ficiently strong motional and/or external electric fields 
are investigated in section Ull Dl first at hand of the hy- 
drogen atom. A possible pathway for their systematic 
experimental preparation is outlined. A substantial part 
of this section analyzes the giant dipole states i.e. res- 
onances of multiply-excited Rydberg atoms in crossed 
fields. Section IIII El contains an application of the gi- 
ant dipole physics to matter-antimatter systems specif- 
ically the positronium atom. As a major outcome it is 
shown that giant dipole states of positronium represent 
extremely long-lived and therefore quasistable states of 
matter-antimatter systems [901 ] ■ 

The fundamental properties of moving atomic systems 
in homogeneous magnetic and/or crossed homogeneous 
electric and magnetic fields are of immediate relevance 
to the dynamics of matter as well as antimatter Rydberg 
gases and plasmas in strong magnetic fields. The cou- 
pling of the center of mass (CM) and internal (electronic) 
motion in a homogeneous B field has consequences for the 
electronic structure of the atomic system. The emergence 
of giant dipole states (GDS) due to the double- well struc- 
ture of the generalized potential is one such consequence 
being a result of the gauge-independent pseudoseparation 
of the CM motion. 



A. The Hamiltonian - symmetries and constants of 
motion 

We consider a finite number of interacting charged par- 
ticles exposed to an external homogeneous and static 
magnetic field employing the non-relativistic Coulomb 



interaction. Relativistic effects induced by the external 
magnetic field become relevant only [U [92|, [93[ for field 
strengths well-above the regime considered here, i.e. par- 
ticularly for field strengths occuring in the photospheres 
of neutron stars {B 3> 10 5 T). The starting-point of our 
analysis is the nonrelativistic Hamiltonian in Cartesian 
coordinates in the laboratory frame 

^ = £^ + £^-^•1) (i) 

i 1 i>j 

where we have omitted the interaction of the spins with 
the magnetic field. The latter is indispensable to obtain 
the right ordering of the energy levels in the astrophysi- 
cally strong field regime 0, [9jt H [97|, [98| but plays no 
role for the following considerations focusing on the spa- 
tial symmetries and constants of motion and the resulting 
transformations to the Hamiltonian. We remark that the 
following analysis holds not only for the case of V c rep- 
resenting the Coulomb potential but for any translation- 
and rotation invariant potential. ~ki is the kinetic (me- 
chanical) momentum for the i-th particle of mass m*, 
charge qi, position Tj with the canonical momentum p^, 
and the vector potential A^ 

7T,; = p s ; - q t Ai (2) 

The Schrodinger equation Hl^l = E*f? l belonging to 
the Hamiltonian in Eq. ([1]) is invariant with respect to 
the gauge transformations 



* L i — ► exp ^+i£ (?i A(r i )J * L (3) 

A< i — ► A,-+ViA(ri) 

In this and the following section we will, as far as 
possible, use a gauge-independent formalism. In those 
cases where it is necessary to introduce an explicite 
gauge, we make use of the symmetric Coulomb gauge, 
Ai = | (B x rj), where B is the magnetic field vector. 
The kinetic mechanical momentum ^ obeys the com- 
mutation relation 

[-K ia ,TVjp] — iq^apjBySij (4) 

where the Greek indices denote the components of the 
corresponding vectors. e Q /3 7 is the completely antisym- 
metric tens or of rank three. The so-called pseudomomen- 
tum ki [ai, [Tool GUI 

ki = TTi + qi [B x n] (5) 

obeys the following commutation relations 

\ki a ,Tjp] = (6) 
[ki Q ,kj£j] = iqi£a/3fBj5ij 
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The mathematic al and ph ysical interpretation of the 
pseudomomentum jlOdl . Il0l| depends on the situation 
under consideration and is closely connected to the phase 
space symmetries and conserved quantities of the under- 
lying Hamiltonian. In the symmetric gauge, the kinetic 
and pseudomomcnta differ only by the sign of their field- 
dependent terms. The classical Hamiltonian correspond- 
ing to the Eq. ([T]) is then invariant under the phase space 
translation 



k ± = qr c x B; r c± = — (Bxk) (12) 

This extends naturally to a system of noninteracting 
particles where the total pseudomomentum is interpreted 
as the superposition vector of the guiding centers of the 
individual particles. For a neutral two-body system, we 
have 



(r i ,p i )^(r i + a, Pi + (|)Bxa) (7) 

which is a combined translation and gauge transforma- 
tion. In the absence of the external magnetic field, the 
phase space translation group reduces to a translation 
in coordinate space forming an Abelian group for both 
underlying classical and quantum systems. The coordi- 
nate translation group is generated by the conserved total 
canonical momentum P = p; being equal to the total 
kinetic momentum. The phase space translation group 
in the presence of the field for a quantum system is non- 
Abelian, i.e. phase space translations do not commute 
in general. The generators of the phase space translation 
group are the components of the total pseudomomentum 

Kl = £> (8) 

i 

which represent constants of motion in the presence of 
the field 



K L = q (r cl - r c2 ) x B (13) 

i.e. Kx is proportional to the distance vector of the 
guiding centers of the two oppositely charged particles. 
The dual interpretation of the pseudomomentum be- 
comes evident if one takes i nto a ccount that in case of the 
interaction with radiation [lO0j |. the sum of the photon 
momentum and the pseudomomentum is conserved. In 
the absence of the homogeneous field, but for interacting 
particles, the pseudomomentum K/, reduces to the con- 
served total canonical momentum P leading to an exact 
separation of the CM motion. For an interacting particle 
system in the presence of the field the interpretation of 
the pseudomomentum possesses an 'interpolating charac- 
ter' between the kinetic and positional representations. 
The guiding center approach to the classical dynamics 
is adequate in the regime where the forces due to the 
magnetic field dominate the Coulomb forces (see section 

Unfit . 

A further ex act consta nt of motion of the Hamiltonian 
is provided by pM OH 



[Hl,K l }=0 (9) 

As indicated above, the components of the pseudomo- 
mentum Ki can, in general, not be made sharp simulta- 
neously, 




In case of the symmetric gauge for the vector potential, 
this quantity is identical to the projection of the total 
orbital angular momentum onto the magnetic field axis 



[K Lct ,K L/3 ] = -iQeapyBy (10) 

where Q = J2i 9» i s the total charge of the system. 
Only in the case of a neutral particle system, i.e. for 
Q = 0, the components of Kj, commute altogether and 
can be exploited equally as constants of motion in the 
context of a pseudoseparation of the CM motion. We 
remark that since [K L \\, Kl±] = 0, the CM motion along 
the B field is exactly decoupled from the relative motion. 

For a single charged particle in a homogeneous mag- 
netic field, the coordinate vector of the center of the clas- 
sical Landau orbit is provided by the so-called guiding 
center, see a discussion of the guiding center approxima- 
tion later in Sec. IIII Fl 

r c = r+-l(^xB) (11) 

The guiding center vector is closely related to the pseu- 
domomentum vector 



£ ii = b b E^ x p^ ( 15 ) 

i 

£|| is the generator for rotation about the magnetic 
field axis, and does not commute with the components 
of the pseudomomentum perpendicular to the magnetic 
field but obeys 

[£||,X L ||] = [£||,KiJ =0 (16) 

As a consequence, a complete set of commuting con- 
stants of motion for a neutral system, would either be 
the components of or {K L ^ , K| ±) £||}. For charged 
systems comprising of heavy and light particles, such 
as atoms and molecules, approximate constants of mo- 
tion exist. These quantities are particularly helpful for 
the transformation of the Hamiltonian such that its de- 
pendencies on the CM and internal degrees of freedom 
take on a particularly simple and distinct appearance 

[ MIMfiMfM lio&l. 
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B. Pseudoseparation of the center of mass motion 

The pseudomomentum is closely related to the center 
of mass motion of the underlying particle system. To 
elucidate this fact, let us perform a transformation from 
the laboratory coordinate system to CM and internal co- 
ordinates. In the following we focus on atomic systems 
consisting of a heavy nucleus with mass Mq and N elec- 
trons each with mass m. For the symmetric gauge, the 
pseudomomentum then takes on the following appear- 
ance 



K = P+^BxR+^aBx^r, (17) 

i 

where (R, P) are the CM coordinate vector and its 
canonically conjugate momentum, respectively. {ri,pi} 
denote the relative coordinate vectors of the electrons 
with respect to the nucleus and their canonically con- 
jugate momenta, respectively, a — (Mo + Zm)/M,Q = 
(N — Z)e and Z, M, e are the nuclear charge number, 
the total atomic mass and the electron charge, respec- 
tively. For Q — 0, i.e. a neutral atom, the pseudomo- 
mentum does not depend on the CM coordinates. In this 
case all components of the pseudomomentum commute 
i.e. they can be made sharp simultaneously. Introducing 
the pseudomomentum as a canonical momentum elimi- 
nates the cyclic CM coordinates from the Hamiltonian 
thereby resultin g in a so-called pse udoseparation of the 
CM motion [MEM EM EM Ell ■ We emphasize that 
the pseudoseparation does not yield a complete separa- 
tion of the CM and internal motion. Indeed, an intricate 
coupling of the CM and internal motion remains, leading 
to effects such as the classical diffusion of the CM for a 
highly excited Rydberg atom (see section Ull CI) . 

To elaborate the case of a neutral atom in more detail, 
we observe that the eigenfunctions of the Hamiltonian 
can simultaneously be chosen as eigenfunctions of the 
pseudomomentum. Since the pseudomomentum is linear 
with respect to the coordinates and momenta in Eq. (|17[) . 
the dependence of the eigenfunctions on the CM has the 
following form 

*({r,},R;/C) = exp^/C-|Bx^r^ rJ tf ({*};£) 



e X p(i/CR) t/ n *o ({r;};/C) 



(18) 



with the eigenvalue equation ICW = KSS> . Apart 
from the plane wave component, exp(i/CR), the unitary 
transformation U n represents the only CM coordinate- 
dependent part of the total wave function The wave 
function \&o of the internal motion depends paramet- 
rically on the eigenvalue K, of the pse udomomen tum. 
Transforming the Hamiltonian with U n |l0ll . Illdl | and 
replacing the pseudomomentum with its eigenvalue, pro- 
vides us with the following pseudoseparated Hamiltonian 



«-5B*"-i(«:xB)V B 



where Hi nt represents the part of the Hamiltonian 
which depends exclusively on the electronic degrees of 
freedom. The first term of Eq. (fT9|) is a constant energy 
shift. The second term represents the coupling of the CM 
to the electronic motion. It is a Stark term [99j involving 
a motional electric field -g/C x B: due to the collective 
motion of the neutral atom in the magnetic field the atom 
feels an additional homogeneous electric field perpendic- 
ular to the magnetic field axis which causes a separation 
of the centers of charge of the atom, i. e. it forms an 
electric dipole. 

There have been a number of investigations addressing 
the effects of the motional Stark field or finite nuclear 
mass effects in general for low-lying electronic states of 



atoms in strong magnetic fields primarily for the two- 
body problem i.e. the hydrogen ato m [l 
EH El EH EM EH flia. fll9L fl20j a nd more recently 



ly fo r the two- 

uM EH EH, 



also for the helium atom |l2ll I122L Il23l Il24j . [For early 
experimental work on the influence of motional Stark 
field on the spectra a nd t r ansit ion line shapes, we refer 
the reader to Refs. jl25l Il26j ]. Furthermore, a num- 
ber of early work relating to various two-body aspects of 
part i cle systems in cros sed electric and magnetic fields 
[1271 EH, EM EM EMI , are worthy of mention. 

The pseudoseparation of the CM motion for neutral 
systems yields a natural division of the Hamiltonian into 
three separate parts involving the CM, its coupling to the 
electronic motion and the electronic degrees of freedom, 
respectively. The transformed Hamiltonian was obtained 
for a specific symmetric gauge. A ma nifestly ga uge- 
independent approach has been developed [l32l . |X33( ] pro- 
viding us with a gauge-independent potential picture for 
moving atoms in magnetic fields. 

Writing the general gauge as A(r) = x r + VA(r), 
the first step is to perform the coordinate transformation 
from laboratory to the CM and internal coordinates. By 
constructing the common eigenfunctions of the Hamil- 
tonian and the pseudomomentum, and using the specific 
dependence of the gauge function A on the new variables, 
it becomes possible to derive a unitary transformation 
that eliminates the CM coordinates from the Hamilto- 
nian. Rearranging terms, we arrive at the Hamiltonian 
H = T + V with 



T 



+ 



l — J2 (pi- | B x r » + e P B x E r * + eV ^) 

\ i i i / 

^(/C-iVaBx(lE-)) 2 



2M 

eE^r, + K({r,}) 



(21) 



Hi, 



(19) 



where (3 = ||,7 = In the above Hamilto- 

nian, we have additionally taken into account an external 
homogeneous electric field E. /({r^}) is an integration 
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constant resulting from the differential equation that en- 
sures the total wave function to be an eigcnfunction of 
K. T represents the kinetic energy of the electrons in 
the presence of the magnetic field and this term is, as 
expected, explicitly gauge-dependent via the scalar func- 
tion/. V, on the other hand, is independent of the chosen 
gauge (no scalar function / occurs in V) and can there- 
fore be interpreted as a generalized potential. Besides 
the Coulomb interaction terms V and the electric Stark 
term due to the external electric field the first quadratic 
term of the potential V is of particular relevance. Apart 
from the trivial constant K? /2M it contains a motional 
electric field term e/M(B x JC) r i an d a diamagnetic 
term e 2 /2M(B x ^ r^) 2 . The relevant quantity occuring 
in the latter two potential terms is the electronic center of 
mass (ECM), i.e. R = A r * m the internal coordinate 
frame. It is therefore the ECM which experiences inter- 
actions beyond the Coulomb potential and which enters 
the generalized potential for multi-electron systems. In 
case of (effective) one-electron systems the ECM reduces 
(approximately) to the coordinate vector of the single 
electron. 

We remark, that the first quadratic term in Eq.lJ^D]) 
which is an important part of the total potential V, rep- 
resents the kinetic energy of the CM of the atom, which 
can be verified by inspecting the Hamiltonian equation 
of motion for the CM 




Therefore the CM kinetic energy of a neutral atom pro- 
vides a potential for the internal motion of the electrons 
of the atom. This kinetic energy is due to the vanish- 
ing net charge of the system independent of any chosen 
gauge of the vector potential. 

For completeness, we give here the Hamiltonian Hi = 
T\ + Vi for the hydrogen atom. It emerges from Eqs. 
(|20|2ip by specializing to N = 1, subsumming the two 
quadratic terms of the kinetic energy in Eq. (f2T)|) , as 
well as introducing the reduced masses /i = and 

/ _ mM 
r M -m 

7i = ^(p-^Bxr + eV/W) 2 

Vi = ^-(/C-eBxr) 2 + V c (r)-eEr (23) 

The generalized potential Vi gives rise to a number of 
intriguing features such as a double well potential which 
accommodates weakly bound Rydberg states with huge 
electric dipole moments, see also Sec. (jlll Dp for a 
discussion of GDS for multi-electron atoms. Here we 
wish to discuss in some more detail the properties of 
the generalize d potential Vi for a single-electron atom 
[132[ Il34l Il35l | . In Figure [H we show a two-dimensional 
intersection of Vi in the [x, y) plane perpendicular to the 
magnetic field. The saddle, due to a competition between 




FIG. 2: A two-dimensional intersection of the potential Vi in 
the (x, y) plane perpendicular to the magnetic field 



the Coulomb potential and the linear Stark term, is ev- 
ident. At larger values of x, the quadratic diamagnetic 
potential begins to dominate, leading to the formation of 
a shallow and long-range potential minimum. Discrete 
energy states in this potential well are the GDS, see also 

Sec. Imp] 

The existence of the saddle point and the outer well, 
depends, on the strength of the B field and the pseudomo- 
mentum or, respectively, on the strength of an external 
electric field oriented perpendicular to the magnetic field. 
Indeed, the Stark term due to the external electric field 
E may be combined with the quadratic term that repre- 
sents the kinetic energy of the CM motion in Eq. ([23]) by 
redefining the pseudomomentum K, — > K, + Mvd where 
v<j is the drift velocity of charged particles in the crossed 
external fields (E, B). The equilibrium position of the 
potential well is obtained from the following conditions: 

y = z = ; x 3 +(K/B)x 2 -(M/B 2 ) = and x < 

(24) 

The condition for the existence of both the saddle point 
and the outer well can be stated as follows: 

27 

K 3 > —BM (25) 
4 

The position ro of the minimum is then obtained from 
Eq.®, 



yo = z Q = : x « (-K/B) + KM / (K 3 - 2MB) 

(26) 

Due to the confining behaviour of the generalized po- 
tential Vi perpendicular to the magnetic field, it is evi- 
dent that ionization of an isolated atom in a giant dipole 
state can only take place parallel to the magnetic field. 
The quantized states in the outer well (GDS) are dis- 
cussed in Sec. IIIID1 
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C. Center of mass coupling to the internal 
dynamics 

The pseudo-separation of the CM motion and adoption 
of the symmetric gauge for the hydrogen Hamiltonian in 
a magnetic field results in 



Hi = 



T1 + V1 
If e fi 



2fi 



2/i' 



Bxr 



(27) 

^(/C- e Bxr) 2 + V c 



- 27^ + C (Bxr)2 + K 



(28) 



The B field is along the z-axis. In the first equation, there 
is a clear division of kinetic and potential energies, the 
term containing K is an effective potential for the inter- 
nal motion, the second form of the Hamiltonian does not 
have the usual division of the kinetic and potential en- 
ergies and contains the linear Zeeman and the quadratic 
diamagnetic terms. The kinetic energy of the CM motion 
is, as already indicated in section pil Bp in the context 
of Eq.(f2"2"]). provided by the quadratic term of the gener- 
alized potential for the internal motion i.e. we have 



f 



f 



— — R = — (K-eB x r) 



2M 



2M 



(29) 



To understand the effects of the coupling between the 
CM and electronic motion we focus in the next two sub- 
sections on the classical dynamics of the hydrogen atom. 
Although being a quantum (or semiclassical) object a 
classical description of the highly excited atom is illumi- 
native and provides relevant insights i.e. a qualitative 
description of many of its properties can be obtained via 
a classical picture (see refs. [EH, [EH [E3, [EH and refs. 
therein) . 



1. The classical center of mass motion for a vanishing 
pseudomomentum 

The Hamiltonian equations of motion for the CM and 
electronic motion for the hydrogen atom are 



K, = 

R = —K- — (B x r) 

M M ' 

1 e .„ 
r = -p- — B x r 

/t 2/r 

p = 4(B^/C)~(Bxp) 



M 



2/i' 



— Bx (Bxr)-e 2 ^ 
4/i |r| 3 



(30) 
(31) 
(32) 

(33) 



The CM velocity couples to the relative coordinates, 
transverse to the B field, even when K = 0. The CM 



equations of motion are independent of the gauge. The 
Hamiltonian equations of motion for the electronic de- 
grees of freedom (|32I33|) are for the case IC = indepen- 
dent of the CM motion, i.e. they decouple from the CM 
motion. 

For IC = 0, S||z-axis the Hamiltonian possess an az- 
imuthal symmetry i.e. the projection of the orbital an- 
gular momentum onto the z— axis is a conserved quan- 
tity. Let us focus here on the case L z — 0. Apart from 
similarity transformations, the classical dynamics does 
not depend on the field strength and energy separately 
but only on the scaled energy e = EB~ 2 / 3 . Varying e 
from —3 to —0.1 the classical phase space belonging to 
the electronic motion of the hydrogen atom undergoes 
a transition from (almost) complete regularity to a fully 
chaotic situation (see Ref. [54[ and refs. therein). 

The natural question arises what kinds of center of 
mass motion are possible if the internal motion goes 
through the whole range from regularity to irregularity. 
Another closely related question results from the follow- 
ing considerations. The internal motion is restricted to 
the energy shell and takes place in a bounded region of 
phase space. The phase space of the CM motion how- 
ever is, at least in principle, unbounded and therefore 
one may ask whether or not phase space is filled out by 
the CM trajectory and how this depends on the regular - 
ity or irregularity of the internal motion |l36l Il37l Il38l | . 
We remind the reader that in the absence of the external 
field the CM motion decouples from the internal motion 
and exhibits a straigth-lined motion. 

Let us begin our discussion with the 'deep' regular re- 
gion, i.e. the region for which the use of low order classi- 
cal perturbation theory with respect to the diamagnetic 
term is valid. The classical trajectories for the internal 
motion of the hydr ogen atom in this region have been 
classified in ref. 139]: There are the so-called rotators and 
librators covering the complete or only part of the inter- 
val [0,27r] for the corresponding angle variable. Figure [3] 
shows a typical CM trajectory if the internal motion is 
of the rotator type. The CM performs an approximately 
smooth circular motion on a ring and is, therefore, lim- 
ited to a small range of coordinate space. According to 
Eq. (f3"Tj) for IC — 0, the quasi-periodic behaviour of the 
electronic degree of freedom r goes hand in hand with 
the quasi-periodic motion of the CM degrees of freedom. 

Let us now turn to the situation of a fully chaotic phase 
space for the internal motion. Figure [4] shows a typical 
CM trajectory for this case. The CM motion is, indepen- 
dently of the special initial conditions, no more restricted 
to some bounded volume of phase space. With increas- 
ing time the CM trajectory fills out an increasing volume 
of the CM hyperplane. Its propagation closely resembles 
the random motion of e.g. a Brownian particle. After 
all, this is quite natural since the underlying equation of 
motion (|3"Tj) for t he CM is a Langevin-type equation with- 
out friction [l40j . The corresponding stochastic Langevin 
force is replaced here by the deterministic chaotic force 
— e (B x r). A characteristic feature of the random Brow- 
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FIG. 3: A center of mass trajectory in the CM hyperplane 
perpendicular to the magnetic field with starting-point R(t = 
0) = (0, 0, 0), for E = -10" 3 , 10~ 5 (B= 2.35 T), L z = 0, e = 
—2.15, all in atomic units. From [1381 ] . 
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FIG. 4: A center of mass trajectory in the CM hyperplane 
perpendicular to the magnetic field with starting-point R(f = 
0) = (0,0,0). The underlying total energy E, magnetic field 
strength B, angular momentum L z and scaled energy e are: 
E = -4.7 • 10~ 5 , 10" 5 (this corresponds to 2.35 Tesla), L z = 
0, e = —0.101. All values in atomic units. From [l38t ] ■ 



nian motion is its diffusion, i.e. the linear dependency of 
the travelled mean square distance on time, 

( P 2 CM ) =D-t (34) 

where Pq M = Xq M + Yq M and D is the correspond- 
ing diffusion constant. The chaotic deterministic mo- 
tion of the CM possesses properties of randomness in 
the sense that it provides the well-known diffusion law 
of random walk models (f or a discussion of randomness 
versus chaoticity see Refs. [HJBIlHISl). Let us briefly 
address the dependence of the diffusion constant on the 
energy and magnetic field strength. With increasing en- 
ergy (always within a completely chaotic phase space) 
the mean square distance of the electron from the mag- 
netic field axis passing through the nucleus increases. As 
a consequence the mean CM velocity and therefore also 
the diffusion constant increases. The diffusion constant 
increases with increasing magnetic field strength thereby 
showing a saturation behaviour. The latter can be un- 
derstood by inspecting the explicit as well as implicit 
field-dependence of the CM velocity given by Eq. (l3"Tj) for 
K = 0. To conclude, the coupling of the CM and internal 
motion due to the external homogeneous magnetic field 
causes a distinct transition in the classical CM motion of 
the atom from bounded oscillatory to unbounded diffu- 
sive motion. We emphasize that this transition happens 
for an isolated atom and is of inherently different ori- 
gin than the random motion taking place for Brownian 
particles which show diffusive behaviour due to random 
collision events among different atoms. 

2. The classical center of mass motion for a nonvanishing 
pseudomomentum 

Let us now study the coupled CM and electronic mo- 
tion for the case of a nonvanishing pseudomomentum. 
Inspecting the equations of motion (|30I31I32I33[) . one ob- 
serves that the CM velocity now contains an additional 
constant term proportional to the pseudomomentum. In- 
tuitively, one would therefore expect a straightlined CM 
motion being superimposed on the motion resulting from 
the corresponding coupling to the internal motion (see 
previous subsection IIII C 1| . Additionally, the equations 
of motion for the internal degrees of freedom depend on 
the value of the pseudomomentum. This is a major dif- 
ference with the case of a vanishing pseudomomentum, 
where the internal motion was decoupled from the CM 
motion. 

A major new consequence for the classical dynamics 
of the highly excited hydrogen atom with nonvanishing 
pseudomomentum is the intermittent behaviour of the 
corresponding trajectories, i.e. the appearance of alter- 
nating phases of quasiregular and chaotic motion [l34| . 
Figure [5] shows for a typical trajectory, the projection 
of the electronic motion on a plane perpendicular to the 
magnetic field axis. One immediate observation is that 
there exists two alternating types of motion. During one 
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phase of the motion, the electron and the nucleus are in 
the x, y— plane close together and this shows up through 
the black dimple close to the origin in Figure [5] During 
this phase of the motion, the Coulomb and diamagnetic 
interactions are comparable and the trajectory is, there- 
fore, chao tic, having a finite local Lyapunov exponent 
[144 , 1 145[ | During the other phase, the electron and the 
nucleus move far apart from each other. The relative 
motion in the x, y— plane then approximately takes place 
on a circle with a large radius. The Coulomb energy 
provides here only a small correction to the strong mag- 
netic interaction. This phase of motion has an essentially 
vanishing local Lyapunov exponent. 

The radius of the circle of the regular phase of the elec- 
tronic x, y— motion of an intermittent trajectory can be 
understood in terms of the pseudomomentum K. Indeed, 
according to Eq. () 13|) , the total pseudomomentum is pro- 
portional to the cross product of the magnetic field vector 
and the distance vector between the two guiding centers 
for free particles. Since the magnetic field strongly dom- 
inates the Coulomb interaction during the quasiregular 
circular motion shown in Figure El we encounter the sit- 
uation of almost free particles and the radius of the circle 
is approximately given by 

r = -d^ |Bx/c| (35) 

i.e. it is completely determined by the magnetic field 
vector and the pseudomomentum. On the other hand, we 
obtain a completely different interpretation of the pseu- 
domomentum if the electron and the nucleus are close 
together. In this case, the Coulomb interaction is com- 
parable to or dominates the magnetic interaction and 
hence the pseudomomentum is approximately the lin- 
ear kinetic momentum of the translational CM motion. 
Fig. [5] shows the CM motion for the trajectory whose 
electronic motion is also given in Fig. It consists 
of alternating phases of purely translational and circu- 
lar motions. As already mentioned, the electron and the 
nucleus are strongly bounded, i.e. close together, during 
the time interval of chaotic electronic motion. This is pre- 
cisely the time period during which the CM performs a 
purely translation motion. The time period of quasiregu- 
lar circular electronic motion corresponds to the period of 
approximately circular CM motion. Intcrmittency there- 
fore shows up in the CM motion by alternating phases of 
more or less straight-lined and circular motion. 

Inspecting the corresponding electronic and nuclear 
motions in the laboratory coordinate system, it turns 
out that the large amplitude motion perpendicular to 
the magnetic field, shown in Fig. [5] in the case of relative 
motion, is performed by the nucleus whereas the electron 
remains essentially localized close to the origin. This pe- 
culiarity has its origin in the fact that the quasiregular 
phases of the intermittent motion are always character- 
ized by large negative L z . Motion with L z « 0, how- 
ever, can, in the laboratory coordinate system, only be 
performed by the nucleus. During the period of large 
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FIG. 5: Typical intermittent trajectory of the electronic mo- 
tion for nonvanishing pseudomomentum: Projection onto the 
x, y— plane - top subfigure, center of mass motion perpen- 
dicular to the magnetic field for R(t = 0) = (0, 0, 0) - 
middle subfigure, and internal z— coordinate as a function 
of time. Field strength, energy and pseudomomentum are 
B = 1CT 5 ,.E = 1.722 ■ 1CT 4 , K = (0,1,0). All values in atomic 
units. 
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amplitude motion of the nucleus perpendicular to the 
magnetic field , the electron performs a large amplitude 
motion parallel to the magnetic field. Fig. [5] shows the 
electronic relative motion parallel to the magnetic field 
as a function of time for the same trajectory whose CM 
and internal motion perpendicular to the magnetic field 
are shown in Fig. [51 respectively. 

When the nucleus returns to the electron performing 
large amplitude motion parallel to the field axis, the two 
remain close together for the next phase of chaotic motion 
close to the origin. During this period of motion the CM 
performs to a good approximation a purely translational 
motion (see FigEJ) . We remark that the large amplitude 
motion perpendicular to the magnetic field circulates, in 
case of Eq. ((25]) , the outer well whereas for the period of 
chaotic motion where the electron and nucleus stay close 
together the dynamics takes place close to the origin i.e. 
in the well created by the Coulomb interaction. 



Hamiltonian reads 7i a = T\ + W. This is the Hamilto- 
nian of a charged anisotropic oscillator in a homogeneous 
magnetic field. Since Eq. (|25f is fulfilled, all frequencies 
belonging to the harmonic potential terms in Eq. (|36[) 
are positive: 

* = («)' (39) 

The spec trum of this oscillator in the B field is har- 
monic |132| 



D. Giant dipole states in crossed fields 

1. The hydrogen atom 

The generalized potential (see Eqs. I21|27|) exhibit, for 
the case of a sufficiently large motional or external elec- 
tric field (see Eq. (j25)) ). a double well structure. The un- 
derlying classical dynamics has been shown above to be 
characterized by an intermittent dynamics. Naturally, 
the question arises whether there exist quantum states 
which live exc lusively in the outer we l l an d what their 
properties are [HI EM EM EM EM EH- To explore 
their properties, let us firstly expand the generalized po- 
tential Vi in Eq. (f2"7) around the minimum in the outer 
well, whose position will be denoted by r = (xo,0, 0) 
transverse to the field which is oriented along the z— axis 
and a motional or external electric field along the x— axis. 
The expansion of V\ around the minimum ro of the outer 
well has to be accompanied by a gauge transformation 
which introduces ro as the new gauge center for the ki- 
netic energy 71 . An alternative equi valent pat h which has 
been outlined and followed in Refs. [lOOl EH is to apply 
a momentum-dependent unitary transformation that re- 
sults in a decentered Coulomb potential. Performing the 
above expansion up to the second order and employing 
the conditions in Eq. (f24j) . we arrive at the following form 
W for our approximation to the generalized potential Vi 
close to its outer minimum 

w = {m + 4) x2+ % -^k) y2 - (4) z2+c 

(36) 

with the constant C = — ^jX^j and (x,y, z) de- 
note now the elongation coordinates from the minimum 
r . W is the potential of a three-dimensional anisotropic 
charged harmonic oscillator. The underlying kinetic en- 
ergy is provided by T\ in Eq. (|27p i.e. our approximate 



E n+n _ riz = Co + yn + + -j uj + + |^n_ + -j 0J- 

+ (n z + ^ju! z (40) 

with Co = (-^ — ym x o 2l7^ 2 ) an< ^ ^he normal 
mode frequencies read 




with the cyclotron frequency ui c = —^f-. The ener- 
getically low-lying CDS in the outer potential well be- 
ing well-approximated by the eigenstates of the charged 
anisotropic harmonic oscillator in a magnetic field are 
therefore inherently different from the Rydberg states in 
the Coulomb well. In particular, they possess a huge 
dipole moment along the electric field axis which is of 
the order of —exo. We remark that for large values of K, 
the ionization threshold is much closer to the energy of 
the minimum of the outer well than to the corresponding 
energy of the saddle point. 

To quantify the deviation from the harmonic spectrum 
with increasing degree of excitation in the outer well Fig [5] 
shows the energy difference of the eigenvalues according 
to Eq. (|4i"j) and the exact eigenenergies of the GDS of the 
hydrogen atom in the outer potential well in units of U- . 
This difference grows stepwise while neighboring states 
show very different deviations from the harmonic approx- 
imation. To explain these features let us look e.g. at the 
energy level 331. We see that the difference between ex- 
act and and approximated energies for this level is much 
larger than for the levels below 331. The level 331 has 
the quantum numbers n + = n_ = and n z — 10, i.e., 
the quantum number n z = 10 appears for the first time. 
Looking at higher levels, there are maxima of energy dif- 
ferences AE/uj- every 11th level above 331 up to level 



13 



3»» 



200 
level 



FIG. 6: Energy difference between the eigenvalues of the 
anisotropic charged harmonic oscillator in a magnetic field 
and the exact eigenenergies of a hydrogen atom in crossed 
electric and magnetic fields in the outer potential well in units 
of ui- as a function of the energy level. Pa rameter values are 
B = 10" 5 a.u. and K = 0.6 a.u. From [l3l. 



397. For these levels (342, 353, etc.), the quantum num- 
ber n z — 10 and n_ = 1,2, ... Between two levels with 
n z = 10 there are levels with n z < 10, and apparently 
the energy difference for these is smaller. Therefore, the 
difference between harmonic and exact energies is mostly 
determined by the quantum number n z . Hence, the an- 
harmonicity of the exact potential is most pronounced in 
the z-direction. This can also be seen in perturbation 
theory for higher terms of the expansion of the Coulomb 
potential where the major contributions to the energy 
corrections are due to those terms containing high pow- 
ers of z. 

In view of the novel properties of the GDS, schemes 
for their experimental realization are also interesting. 
There have been experiments that indicate the existence 
of atoms with very large dipole m oments in crossed elec- 
tric and magnetic fields |l49l . 1 1 50t | for energies above the 
saddle-point energy. These were performed in compara- 
tively weak fields, B = 10~ 6 a.u., and the resulting esti- 
mate for the dipole moment of the Rydberg atoms was 
1.9 x 10 4 Debye. This is roughly what one should expect 
for the value of the dipole moment for saddle point states. 
The experimental technique employed to determine the 
atomic dipole moment was to add a slight inhomogene- 
ity to the electric field and to measure the deflection of 
the atoms, this deflection being proportional to the field 
inhomogeneity and the atomic dipole moment. 

Finding deeply-bound GDS is a large experimental 
challenge. Starting from the ground stat e, th is requires 



a transition via Rydberg states to GDS 151]. Alterna- 
tively one could think of starting with free particles in the 
continuum and forming GDS via three-body recombina- 
tion in crossed fields, which is a typical process occuring 
in ultracold dense plasmas (see section IPTf . The key idea 



to 'transport' the quantum states of the atom from the 
centered ground to the decentered GDS is to combine 
a laser excitation to a Rydberg state in the Coulomb 
well, with a two-step switching procedure of an exter- 
nal electric field which adiabatically transfers the quan- 
tum states to a narrow distribution of strongly bound 
GDS. This scheme is based on the switching of an ex- 
ternal electric field whereas the corresponding motional 
electric field should be negligible, as is the case e.g. in 
an ultracold environment. Following the laser excitation 
to a Rydberg state in the Coulomb well this switching 
procedure is provided by the time-dependent potential 
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The value of the intermediate electric field E c = 
10 _8 a.u. » 5656F/m is reached after t\ = 2.3 x 10 8 a.u. w 
5.6ns. The final value of the external electric field 
is Ef = 2.5 x 10~ 8 a.u. « 12855F/m. This value is 
reached at tf = 9.2 x 10 9 a.u. « 222ns after a delay of 
t 2 - h = 5.7 x 10 8 a.u. w 13.8ns. The period t 2 - h 
of a constant field E c as well as the subsequent second 
switching process ending with the final value Ef of the 
electric field are chosen to achieve as narrow as possible 
a final state distribution in the outer potential well. A 
single switch to achieve the transfer from the Rydberg to 
the GDS, would end up in a broad final distribution of 
the populated states. 

The initial (t = 0) potential consists exclusively of a 
diamagnetic and Coulomb potential terms (see Eq. (|4"2")l 
and Fig. [7]). To simulate the dynamics in the course of 
the switching process we choose an ensemble of trajec- 
tories according to a single chaotic trajectory for a typ- 
ical energy of a Rydberg state in the (quantum-) chaotic 
regime thereby ensuring a random sampling of the energy 
shell. The first switching-on process from zero to the elec- 
tric field strength E c within the time interval < t < t\ 
leads to the formation of a Stark saddle point and a very 
shallow outer well, as Fig. [71 

One needs to keep the electric field constant for a pe- 
riod of time long enough to enable most of the atoms to 
cross the Stark saddle once. The positions of trajectories 
in the (x, y) plane just before the second switching on are 
shown in Fig. [H] A substantial part of the trajectories has 
passed to the outer well by that time. They are seen as 
points scattered around an outer well trajectory. Only 
a tiny fraction of trajectories lead to ionization in the 
course of the switching on process. Then one can increase 
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FIG. 7: Potential energy of the relative motion along the 
external electric field direction. The upper curve corresponds 
to zero electric field [E(t — 0) = 0]. The curve with a shallow 
well corresponds to E = E c (t\ < t < £2). The curve with 
a deeper well corresponds to E(t = tf) = Ef. The mean 
energy of the trajectory ensemble at times t = t% and t = tf 
is shown by dashed-dotted lines joining the turning points 
of the appropriate potential energy curves. The ioni zatio n 
threshold at E = Ef is shown by a dashed line. From [l5lj . 



the electric field up to its final value (see FigJTJ) . In such a 
procedure, the value of E c determines the most probable 
final energy in the outer well. The second switching on 
has to be sufficiently slow in order to reduce the broad- 
ening of the energy distribution during the trapping of 
the relative motion in the outer well. Quantum mechan- 
ically speaking, such a broadening occurs due to transi- 
tions between the adiabatic outer well states. It cannot 
be fully avoided since in the course of trapping the rel- 
ative motion crosses the separatri x be tween the energies 
above and below the Stark saddle [152| . As shown in Fig. 
IH1 all nonionizing trajectories have been brought to the 
outer well at time t — tf. Whether the finally result- 
ing trajectories are truely bound ones in the outer well 
depends sensitively on the parameters chosen. Fig. ITD1 
shows the probability distribution as a function of energy 
for the final ensemble for B = 6 x 10 _5 a.u. It represents 
a distribution with a single dominating peak and a sub- 
stantial portion of the energy being below the ionization 
threshold. The peak is due to the trajectories which had 
crossed the Stark saddle before the second switching on 
phase of the electric field began. The trajectories that 
passed to the outer well during the action of the second 
field switch show up in the energy distribution as a long 
tail for large energies. 



FIG. 8: Trajectories projected on the (x, y) plane at time 
t = t 2 . From QUI. 
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FIG. 9: Trajectories projected on the (x, y) plane at time 
t = t f . From [l5l| ]. 



2. Multi-electron atoms 

We start by discussing the generalized potential in 
Eq. (|2Tj) for a multi-electron atom. Apart from the 
trivial constant K? /2M, it gives rise to a motional elec- 
tric field term e/M(B x JC) r i an d a diamagnetic term 
e 2 /2M(B x J2 r i) 2 - The relevant quantity occuring in 
the latter two potential terms is the electronic center of 
mass (ECM), i.e. R e = jj r » m * ne internal coordinate 
frame. It is therefore the ECM which experiences inter- 
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FIG. 10: Probability distribution as a function of energy in 
the outer well at t > tf for B — 6 x 10~ 5 . The ionization 
threshold is shown by a vertical dashed line. From |l5l| . 



actions beyond the Coulomb potential and which enters 
the generalized potential for multi-electron systems. In 
case of one-electron sy stem s the above potential reduces 
to the one derived in 132] and in particular the ECM 
reduces to the coordinate vector of the single electron. 

Let us focus on doubly excited two-electron systems. 
Since major parts of the generalized potential V depend 
only on the ECM and since both electrons are assumed 
to be highly excited it is natural to introduce the ECM as 
a new coordinate vector. Additionally, we require that 
the kinetic energy should become as simple as possible 
which leads to the relative vector of the two electrons as 
a good choice for the second coordinate vector, i.e. R e = 
(ri + r 2 )/2; r = ri — r 2 . The transformed Hamiltonian 
Ti.2 = T 2 + V-2 therefore decomposes into 
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As can be seen from 



M -2m ' 

the coordinate transformation decoupled the 
kinetic energy terms belonging to the two electronic vec- 
tors and simplified the field-dependent potential terms in 
Eq. @2). 

Doubly excited configurations corresponding to the 
extrema of the six-dimensional potential V 2 (R, r) in 
Eq. (|4"3")) are good candidates for resonances. The roots 
of the six nonlinear coupled equations d V 2 /dr = and 



<9V 2 /<9R = are therefore of immediate interest. With- 
out loss of generality, we assume in the following again 
that the magnetic and electric field vectors point along 
the positive z— axis and negative a;— axis, respectively. As 
previously indicated the Stark term due to the external 
electric field can be taken into account by redefining the 
value of the pseudomomentum. This yields the following 
geometrical conditions 
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where R e = (X e ,Y e , Z e ), r = (x, y, z) and r — |r|,_R e = 
| R e |- Accordingly, r R e = 0, the ECM and the in- 
terelectronic coordinate vector are orthogonal. Since 
Y e = Z e = 0, this leads to x = 0. Furthermore, the 
condition R e — |v3r and equally |R e — ^r| = |R e + ir| 
leads to the fact that the two electrons and the nucleus 
form an equilateral triangle. The remaining nonzero co- 
ordinate X e has to fulfill the corresponding polynomial 
equation P(X e ) = in Eg . (|4"3"| . This completes the 
specification of the extremal configurations which are lo- 
cated on a one-dimensional circular manifold. The elec- 
trons form a decentered triangular configuration and are 
highly correlated through the fact that they are forced 
to stay on opposite sides of a circle. The geometry of 
the extremal configuration described by the above con- 
ditions is illustrated in Fig. [TT] in which the circular ex- 
tremal line as well as the opposite electrons are indicated. 
Both electrons are for laboratory field strengths located 
far from the nucleus, the electron-nucleus distance scales 
with oc -g. The position of the extrema will in the fol- 
lowing be denoted by ro, R e o- 

If the inequality K 3 > |\/3MB is fulfilled, P(X) = 
has two real solutions on the negative x— axis (the de- 
centring direction of the atom) . The smaller of these two 
values (excluding the sign) corresponds to a maximum 
(saddle) of the intersection of the potential V 2 along the 
X— direction, whereas the larger value yields an outer 
minimum. Of particular interest is, of course, the case 
where the ECM is captured in the outer minimum. To 
investigate stability criteria for the case when the ECM 
is located in the outer well, a normal mode analysis is 
necessary in the presence of the external field. The nor- 
mal mode analysis using the second order expansion of V 2 
around the extremal positions yields the following eigen- 
value problem for the harmonic frequencies flj, as the 
energies of the doubly-excited resonances in crossed fields 
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FIG. 11: Shown is a sketch of the geometrical configuration 
of the giant dipole two-electron resonances. The electric and 
magnetic field vector point along the negative x- and positive 
z-direction, respectively. R e is the electronic center of mass 
coordinate and r the relative coordinate vector of the elec- 
trons. The big circular loop lying in the yz— plane indicates 
the g eometrical position of all extremal configurations. From 
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FIG. 12: Eigenmodes {7 P } P for N — 2 electrons as a function 
of the electric field E = BK/M over the range K/K cr £ 
[1, 10] (all modes are imaginary). The three sets of curves 
refer to B = 10" 5 , 1CT 4 , 10" 3 a.u. (from left to right). Each 
top horizontal line represents two degenerate cyclotron modes. 
Below, the two Coulomb modes fall off quickly and intersect 
the CM mode (the nearly horizo ntal line about four orders 
below the cyclotron modes). From |l33l |. 
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and Vj = (Vi, ...,Ve)i are the six-dimensional eigen- 
vectors. The reader notes that the matrix on the r.h.s. 
of Eq. lHH) depends explicitly on Oj which is due to 
the appearance of T 2 in Eq. (|43)l in the presence of the 
magnetic field. 

6 

The frequency spectrum reads as E = ^M-^i + 1) + 

V2(ro,Ro) An analysis of fli finds that the two largest 
frequencies are almost degenerate and are of the order of 
half the electronic cyclotron frequency oj c . We call these 
modes in the following cyclotron modes. Two of the re- 
maining three frequencies belong to the motion parallel 
to the magnetic field which is governed exclusively by 
the Coulomb interaction, the so-called Coulomb modes. 
The remaining frequency is due to the heavy particle dy- 
namics i.e. the CM mode. The frequencies fii are dif- 
ferent by several orders of magnitude and moreover they 
are real for typical laboratory field strengths. We there- 
fore encounter no decay within our harmonic analysis of 
the corresponding resonances which indicates that they 



should possess a significant life time. 

Fig[l2] shows the dependencies of the five nonzero fre- 
quencies on both the electric as well as the magnetic field 
strengths. With increasing magnetic field strength all fre- 
quencies increase. For those frequencies which are associ- 
ated with the cyclotron and CM modes this behaviour is 
evident. For the frequencies associated with the Coulomb 
modes parallel to the magnetic field it is a consequence 
of the fact that the position |Xrj| of the outer minimum 
decreases strongly with increasing field strength. The 
Coulomb potential becomes then stronger and the fre- 
quency in the corresponding well parallel to the field in- 
creases. As can be seen in Fig|T2]the dependence of the 
frequencies on the electric field strength is twofold: The 
frequencies associated with the cyclotron and CM mode 
show in general only a very weak dependence on the elec- 
tric field strength whereas the frequencies associated with 
the Coulomb modes generically exhibit a strong depen- 
dence on the electric field strength (see Fig. [T2|) . 

The above normal mode analysis of the giant dipole 
resonances provides evidence that the two-electron case 
is locally stable, apart from the singular horizontal con- 
figuration. Beyond this analysis, recently a numerical 
ab i n itio study of the resonances has been performed 
[l53l Il54| . Such a six-dimensional resonance study is 
both methodologically and computationally demanding, 
especially in view of the fact that our system is gov- 
erned by dramatically different time scales: The antic- 
ipated order of magnitude is a few picoseconds (ps) for 
the electronic cyclotron motion in strong magnetic fields 
(B w 10~ 4 a.u.) and a few thousand up to a few ten 
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FIG. 13: The motion in X in the case K/K cr = 1.1: Snap- 
shots of the motion of the one-particle density px, reflecting 
oscil lations of both center and shape of the wave packet. From 
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thousands ps for the motion parallel to the magnetic 
field (the typical time scale for the CM mode is a few 
thousand ps). To meet the requirements the multicon- 
figuration time-de pendent Hartree (MCTDH) method 
[I55lll56l . ll57lfl58| has been employed, known for its out- 
standing efficiency in high-dimensional applications. The 
MCTDH approach is a multifunctional multidimensional 
wave packet propagation method which allows to inves- 
tigate both time-independent as well as time-dependent 
problems. Applying it to a fermionic system like ours 
finds its sole justification in the fact that the spatial sep- 
aration between the electrons is so large that they are 
virtually distinguishable. 

Both the stability as well as the spectral properties 
of the giant dipole two-electron resonances have been 
studied this way. Moreover, following the evolution 
of wave packets with an initial displacement from the 
extrem um, t he r obustness of the resonances has been 
tested [HI [HH- For K = l.lK cr , i.e. very close 
to the critical point of the existence of the outer well 
(K cr — ^v3MB), the observed motion display an in- 
stability with respect to some of the degre es of freedom. 
The original wave packet (see Ref. jl53l . Il54| for its prepa- 
ration) shows minor oscillations with respect to the de- 
grees of freedom X, Y which are essentially due to an 
oscillation of its center and shape (see FigJT3]). The rela- 
tive stability in these two directions is readily explained 
in terms of the generalized potential. Apart from a sta- 
bilization by the B field, they experience an additional 
confinement via the quadratic term of the potential V2 
in Eq.(|13]). However, it is evident that the packet broad- 
ens over time, thus slowly delocalizing. Even though the 
relative motion (x,y) perpendicular the magnetic field is 
also gyrationally stabilized, there is no confining term for 
it in the generalized potential V2 as for the ECM. The 



y direction is, for the vertical configuration considered 
here, the zero mode, while the cut through the x direc- 
tion refers to a maximum. In this light, we cannot expect 
the system to be arbitrarily stable in these two degrees 
of freedom. This is a general fact, but for K = l.lK cr 
it is very pronounced. The reason being, near the ex- 
tremum, the Coulomb interaction (the source of these 
instabilities) becomes negligible for higher K. For the 
case K = l.lK cr the zero mode spreading is substantial 
on a time scale of a few thousand ps. The reduced density 
for Z broadens and continually leaks outward. Although 
in the Z direction, the particle lives in a fairly harmonic 
well it is not only lifted in energy to way above the bot- 
tom of the well, but also affected by the instability in x 
via coupling. Turning to the last degree of freedom, z, 
an overall stable behavior is obtained. 

Considering the cases of larger if -values specifically 
K = 2K cr , 10K cr the vertical configuration is found to 
be stable on a time scale of 10ns and 100ns, respectively, 
whereby the latter value was the maximum propagation 
time for K = lOK cr . A closer look at the response of 
the system for K = 2K cr upon displacing Z and z by 
2.000 a.u. unveils that for the excited degrees themselves, 
Z and z, the wave packet is simply reflected between two 
positions ±(z)o with minor (Z) or more pronounced (z) 
deformations and smearing-out due to competing modes. 
The most interesting question may be the effect on the 
{xy) modes. In fact, x is rendered slightly unstable by 
the excitation of the parallel motion. On the time scale 
of several ten thousands of ps the wave packet slowly but 
inevitably starts leaking to outer regions. For the case 
K = 10K cr it is only the y— motion associated with the 
zero mode which shows a minute broadening being the 
first sign of a possible decay which is extrapolated to take 
place for time scale of many microseconds. 

Thus far, we have focused on the vertical two-electron 
configuration, which is defined by the alignment of the 
intcrclectronic vector with the z— axis. Let us now briefly 
comment on the impact of rotations of the two-electron 
configuration around the x— axis. The horizontal con- 
figuration, which corresponds to an alignment along the 
y— axis indeed adds instability, which is discernible even 
for very high K, if less distinct. The in-between diagonal 
configuration is partly unstable on a timescale compara- 
ble to that of the vertical case. Finally, for more than 
two-electron atoms, t he reader is encouraged to follow 
the discussion in Ref. [l53l Il54| . 



E. Matter-antimatter systems in crossed fields 

The lowest order decay rate for positronium T = avp 
in field- free space, where a is the plane- wave cross sec- 
tion for free pair annihilation, v is the relative veloc- 
ity of the electron and positron, and p is the square 
of the wavefunction evaluated at contact, amounts to a 
tenth of a nanosecond (parapositr onium) o r roughly 140 
nanoseconds (orthopositronium) jl59l Il60j |. In view of 
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the transition from a single to a double well structure 
of the potential V and Vi (see Eqs. (|21I23[) .respectively) 
with increasing (motional and/or external) electric field 
strength and the associated emergence of giant dipolc 
states, it is an intriguing question to investigate scenar- 
ios for a combined decentered matter-antimatter system 
[90l ll6lL Il62l . Il63| | . Does the intuitive picture of putting 
one particle into the outer well and the other one to the 
Coulomb well lead to an enhanced stability of these mi- 
croscopic matter-antimatter systems? This seems to be a 
natural conjecture particularly in view of the separation 
of the two wells by a wide and high potential barrier for 
sufficiently strong fields. 

Let us begin by inspecting the underlying Hamiltonian 
for a positronium atom e + /e~ in crossed electric and 
magnetic fields 
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We observe that because for Ps, \i = t^j an d A*' = °°> 
the kinetic energy becomes the free-particle kinetic en- 
ergy, see Eq. (|20|23|) . Second, we encounter a general- 
ized potential term 3^r(K + eB x r) 2 which scales ac- 
cording to oc ^ . This is in contrast to the hydrogen case 
(see Eq. (|23| ). where the scaling is oc jj. However, be- 
cause the minimum of the CM kinetic energy i.e. of the 
generalized potential excluding the Coulombic potential 
is independent of the total mass and the condition for 
the existence of the outer well K 3 > =^-Bm is indica- 
tive that the formation of the well takes place for much 
smaller values of the pseudomomcntum and/or electric 
field strength for Ps atom, much smaller sizes of the GDS 
are possible. Typically GDS of positronium can be of the 
size of a few thousand A (see Fig |T4]) whereas GDS of hy- 
drogen amount to many ten thousands A. 

In Fig.rrj] we show \x \ and \x s \ for different values of 
the B field. The position of the minimum, Xq, increases 
and the position x s of the saddle point decreases with 
increasing value of the pseudomomentum (or the elec- 
tric field strength), common to both the hydrogen and 
positronium atoms. At K = 0.4, B = 5 x 10~ 5 a.u. the 
ionization energy is much closer to the bottom of the well 
than to the saddle point. The barrier width at the ion- 
ization energy is over 7000 a.u. This suggests that the 
bound state spectrum of the exact Hamiltonian should 
be nearly identical to the spectra from the two isolated 
potential wells: as set of Coulomb states localized at the 
origin and a set of outer well states centered around the 
minimum. 

The validity of this isolated potential approach de- 
pends on the amount of tunneling through the barrier, 
or in other words, the amount of mixing between the 
Coulomb states and the outer well states. If the barrier 
is sufficiently high and wide, the eigenstates of the the 
two isolated wells will not mix. The probability density 
will reach nearly zero near the middle of the barrier. If 
positronium is in an outer well state, the electron (or 
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FIG. 14: The location of the saddle point x 3 and the outer 
well minimum x as a function of the pseudomomentum K for 
various values of the magnetic field strength B. The lower side 
of each curve is the value of x s ; x a (K — > 00) — > 0. The upper 
side of each curve is the value of x ; x (K — > 00) — > —00. 
The two segments of the curve coalesce at the critical point 
x c = Xs(K c) = Xo(Kc). All quantities are in atomic units. 
From 1611 



positron) cannot tunnel through the barrier and annihi- 
late. An estimate for the tunneling probability from the 
outer well ground state into the Coulomb well, can be 
obtained from 
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exp 



|p(a;)| dx 



(46) 



where x\ and xi are the classical turning points and 
\p{x)\ = yV(x,0,0) - (wx/2 + C). Lifetimes for tun- 
neling are expected to be of the order of some charac- 
teristic time r « 2tt/lo x divided by the probability P. 
For B = 5 x 10~ 5 and .0878 > K > .0908, the har- 
monic ground state lies above V s ; there are no outer well 
states, only Coulomb states and saddle states which have 
probability in both wells. (Saddle states can occur any 
time the saddle point lies below the ionization energy; 
for certain choices of B and K, all three types of states 
- Coulomb, outer well, and saddle - can be found in the 
spectrum.) The tunneling probability drops dramatically 
with increasing K and the outer well state(s) of positro- 
nium can be considered stable for K > 0.1. The cor- 
responding tunneling lifetime is greater than one year. 
Finite clement calculations fully confirm this simple esti- 
mate of the tunneling probabilities. Giant dipole states 
of e.g. positronium offer therefore a unique pathway of 
preventing microscopic matter antimatter systems from 
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annihilation i.e. for their preparation in long-lived quasi 
stable states. 



F. The guiding center approximation 

In the limit of large magnetic fields the previous analy- 
sis can be greatly simp lified with in the so-called guiding 
center approximation [l64l . 165]. The premise for this 
approximation rests on the decoupling of the different 
degrees of freedom in a strong B field due to their widely 
separated frequencies. If the field is sufficiently strong, 
the cyclotron frequency uj c = eB/m is the largest dynam- 
ical frequency and the cyclotron radius r c is the smallest 
length scale in the system. In this limit, the rapid cy- 
clotron motion can be separated from the remaining de- 
grees of freedom, which then describe the motion of the 
guiding center (see Eq. lfTTjO of the electron. 

Fig. [15] illustrates typical classical trajectories of 
guiding center atoms. The guiding center approxima- 
tion holds when the electron's cyclotron radius is much 
smaller than the distance to the ion. This condition can 
be fulfilled at sufficiently weak binding and also holds for 
the giant dipole states (see section UlI Dp . characterized 
by strongly de-centered electron trajectories in the outer 
potential well. Such trajectories, however, are loosely 
bound to the parent ion, and are, thus, expected to have 
limited lifetimes in collisional environments such as cold 
plasmas. More deeply bound trajectories are confined 
to the central Coulomb well and follow regular orbits 
around the ion. A sketch of the particularly illustrative 
and simp le ca se of a circular, nearly planar guiding cen- 
ter atom |166( | is shown on the right in Fig. 1151 Confined 
by the ionic Coulomb field, the electron's guiding center 
oscillates along the direction of the magnetic field and 
performs an E x B drift motion around the ion. If the 
amplitude z m of the longitudinal oscillations is smaller 
than the radius p of the drift motion around the ion, 
the frequency of the field-aligned oscillations is approxi- 
mately uj z = yj e 2 / (mp 3 ), and the frequency of the E x B 
drift motion is given by u>d — e /(Bp 3 ). Requiring the 
cyclotron frequency to be the largest in the system also 
implies that to z 3> ojd an d leads to the following hierar- 
chy 



u c » w 2 » U>D 



(47) 



These inequalities are fulfilled for large radii of the drift 
motion, p ^> (m/B 2 ) 1 / 3 . At the same time, the hierarchy 
(|47jl implies that the cyclotron radius is much smaller 
than the drift radius p, which is fulfilled if the cyclotron 
energy E cyc is les s than or on the order of the potential 
energy e 2 / p jl66l |. 



1. Equations of motion for guiding center atoms 

To simplify the discussion, let us consider the limiting 
case of a stationary atom M/m — > oo. Within the guid- 




FIG. 15: Left: Classical trajectories of different bound elec- 
tron trajectories for a finite pseudomomentum of K± = la.u.. 
Cyclotron radii have been increased and the c yclo tron fre- 
quency decreased to make them visible. From |l67l ]. Right: 
Schematics of a circular guiding center atom. In descending 
order of the frequencies, the electron performs cyclotron os- 
cillations, oscillates axially along the field and executes E x B 
drift motion around the ion. From Il66ll. 



ing center approximation, the electronic cyclotron energy 
E cyc appearing in the total Hamiltonian 
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represents an approximate constant of motion. Conse- 
quently, the fast cyclotron motion of the electron can be 
averaged out, which yields the following dynamical equa- 
tions for the electron's guiding center momentum and 
position 



mv z = —e 
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y = —bzt 



(49) 



The circular guiding-center atom solution is readily ob- 
tained from these equations, and allows to read off the 
different motional frequencies introduced above. In this 
case r ~ p ~ const, such that the (x, y) motion fol- 
lows a circular orbit with radius p around the ion and 
with the drift frequency wd = e /{Bp 3 )- The frequency 
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uj z = \J e 2 /(mp 3 ) can be directly read off the last equa- 
tion. To leading order in z/p <C 1, the z-motion is de- 
coupled from the remaining degrees of freedom, such that 
one can define the binding energy 



Eh — H z — 



Pz 

2m 



(p 2 + z 2 ) 



3/2 



(50) 



of a GCA atom. If E\> < 0, the electron is bound to the 
ion, irrespective of its cyclotron energy E cyc , i.e. even if 
the total energy is positive. The guiding center approxi- 
mation has proven valuable in gaining useful insights into 
the properties of strongly magnetized Rydberg atoms, 
permitting for exa mple the explicit derivatio n of their 
magnetic moments [l68| . r adia t ive d ecay rates [l69j and 
response to electric fields (1671 Il70l | . The latter will be 
discussed in the next section. 
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FIG. 16: Scaled critical ionization field / = F/(e/p 2 ) of a 
cricular guiding center atom as a fu nction of the scaled axial 
bounce amplitude z m /p. From [1661 ] 
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2. Electric field effects on guiding center atoms 



In sections IIII Dl and MI El some effects of additional 
electric fields on the behavior of strongly magnetized Ry- 
dberg atoms and their utility in preparing and stabiliz- 
ing particular Rydberg states have been discussed. From 
a different perspective, ionization of Rydberg atoms by 
electric fields has proven to be a powerful tool to probe 
their internal states. The magnetic field free case has 
been studied in great detail theoretically and experimen- 
tally (see [l|). Additional magnetic fields can signif- 
icantly complicat e the ioni zation process, and produce 
chaotic dynamics [l7l[ Il72j . On the other hand, in very 
strong magnetic fields the electron motion regains com- 
parably simple character, such that the guiding center 
approximation provides an intuitive picture for field ion- 
ization of strongly magnetized Rydberg atoms. 

Let us consider a circular Rydberg atom with binding 
energy E\> (see Eq. ((50)) ) and M — > oo. The electronic 
guiding center follows a circular drift orbit with radius p, 
and longitudinally oscillates with an amplitude z m <C p, 
determined by 



E h 



Vp 2 



(51) 



The application of a parallel electric field lowers the po- 
tential well that confines the longitudinal z-motion. If 
the field exceeds a critical value F, the electron is able 
to escape and the atom ionizes. For z m = 0, ionization 
occurs as the well vanishes, a nd th e ionization field F is 
straightforwardly obtained as jl67| 



F = a—, 
P 2 



(52) 



For finite axial energies, i.e. for z™ ^ 0, the calculation 
is somewhat more complicated jl66| and has to be per- 
formed numerically. However, one can use some scaling- 
arguments for the axial Hamiltonian H z , and write the 



F 



r/(Wp) 



(53) 



Apparently, the critical field depends only on the drift 
radius p and the relative axial amplitude z m / p, with 
/(0) = a/(4/27). Kuzmin et al. pj56| determined the 
functional dependence of / on the relative amplitude 
from a saddle-point analysis. The result is shown in 
Figdl Eq. ([55|) provides a useful expression to analyze 
state-selective field ionization measurements and will be 
applied in section IV Bl to calculate field spectra of Ryd- 
berg atoms, which were observed in antihydrogen exper- 
iments [zi,[8i|. 

Small electric fields, that do not cause ionization, nev- 
ertheless polarize the atom and thereby affect its center- 
of-mass motion. In situations, where magnetized Ryd- 
berg atoms are formed through recombination in cold 
plasmas (see section [Vj, the plasma's space-charge elec- 
tric field (Fig. ITTa ) c an signific antly alter their motion 
and, as suggested in jl70l . Il73l |. even lead to trapping 
of weakly bound, magnetized atoms. Averaging out the 
rapid cyclotron motion of the electron, one can derive a 
simple and instructive equ ation of m otion for the atomic 
center-of-mass velocity V jl70l . Il73| 



MV = V 



aF 2 



a {Vx • F) (V x B) , 



(54) 



which is valid in the guiding center limit. The electric po- 
larizability a is a function of the magnetic field strength 
B, the drift radius p, the relative axial bounce amplitude 
z m /p and the strength F of the electric field. For F — ► 
and z m <C p , the zero-field polarizability «l is obtained 
as fl7a.ll7a| 
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Fig. [15] compares this simplified approach to the result 
of exact classical calculations that follow the coupled ion 
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FIG. 18: Trapping of a guiding center atom with a initial 
kinetic energy of 4K and binding energy of 40K in the radial 
electric field shown in Fig. 117a . The dashed line shows the 
result of an exact propagation compared to t he p rediction of 
Eq.flSD for a = a(25V/cm) = 1.18a L . From \V7% . 



FIG. 17: Left: Space charge electric field produced by a long 
plasma column of radius 0.2cm and density 2.5 x 10 8 cm -3 . 
Right: Electric field dependence of the polarizability, scaled 
by its zero-field value ql for p = 0.4/im. From [lT3t ] . 

and electron motion and resolve the fast cyclotron os- 
cillations. The figure shows a trajectory of a circular 
atom in an axial magnetic field of B = 3T and the radial 
electric field of Fig. [T7b . as produced by a long plasma 
column. The result of the simplified equations of motion 
(|54")l nicely agrees with the exactly calculated center-of- 
mass dynamics over several orbital periods. Moreover, 
Fig. [TBI shows that already moderate electric fields can 
promote radially bound atom motion for rather large 
binding energies of ~ 40K and center-of-mass energy of 
4K — typical values realized in antihydrogen experiments 
& 

However, such trapping is efficient only for a limited 
range of Rydberg states and kinetic energies [l73j | and 
lacks axial confinement. Several approaches to achieve re- 
liable three-dim ensional confinement have been proposed 
[l74l Il75l ll 76t | . As will be discussed in the next section, 
strongly magnetized Rydberg atoms may, apparently, be 
most effectively confined in magnetic traps, i.e. through 
the application of inhomogenous magnetic fields. 



IV. ULTRACOLD RYDBERG ATOMS IN 
INHOMOGENEOUS FIELDS 

Processing ultracold atoms with static magnetic fields 
requires inhomogeneous field configurations, i.e. spatially 
varying magnetic fields. Since the hyperfine interaction 
is typically much stronger than the interaction of the 



atoms with the external magnetic field, the total angu- 
lar momentum and its resulting magnetic moment pro- 
vide the key quantity which couples to the local mag- 
netic field. Atoms in their electronic ground state and 
some hyperfine substate are within this picture treated 
as point-like magnetic moment carrying particles. For 
most investigations and applications on ultracold atomic 
ensembles in magnetic traps, an adiabatic approach is 
employed i.e. it is assumed that the atomic magnetic 
moment is always aligned with the local magnetic field 
and therefore we can replace the vector coupling /xj-B by 
the scalar one |/xj||B|. This is a simplification which be- 
comes increasingly valid the more the ultracold motion 
of the atoms becomes adiabatic i.e. the change of the 
magnetic field strength seen by the slowly moving atoms 
on the time scale of the precision of their magnetic mo- 
mentum around the local magnetic field should be small 
compared to the local field strength. 

In spite of the frequent use of the adiabatic magnetic 
coupling, theoretical studies of the quantum properties 
of neutral ground state atoms in magnetic traps using 
vector coupling, are of relevance not only for obtain- 
ing the lifetimes but also to reveal potential new classes 
of (trapped or untrapped) eigenstates. Such studies, 
which we will review in the following, have been per- 
formed since the late eighties. Several magnetic field 
configurations have been in the focu s of numerous pub- 
lications, e.^the quadrupole field [177| . the wire trap 
[1781 1 1791 . IX80j ] or the magne tic guide and the Ioffe- 
Pritchard trap [l8l[ Il82l Il83l |. However, only a few 
field configurations allow for stationary solutions. Bound 
states of spin- ^-pa rticles tra pped by a wire have been 
found analytically jl78l . Il84| . The latter authors em- 
ployed a super-symmetric approach in order to derive 
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the Rydberg series of bound states around an infinitely 
thin wire. Finite wire sizes were accounted for by in- 
troducing a quantum defect in the Rydberg formula. A 
numerical derivation of a similar serie s of energies for 
particles with spin up to 2 was given jl8C| by combin- 
ing the finite element method with multichannel quan- 
tum defect theory. However, in most cases atoms are 
trapped in meta-s table states which have finite lifetimes. 



In an early work [1771 ]. dozens of spin-^ resonances in- 
side a three-dimensional magnetic quadrupole field have 
been studied by determining the phas e-shift of scattered 
waves. A similar approach jl8ll . Il82| was chosen while 
studying the dynamics of neutral particles in a magnetic 
guide and corresponding results on resonant st ates were 
reproduced by a complex scaling calculation |l83j . A 
comprehensive study of the spectral properties and life- 
times of neutral spin- g-fe rmio ns in a magnetic guide has 



been performed recently [1851 ]. A wealth of unitary and 
antiunitary symmetries has been revealed giving in par- 
ticular rise to a twofold degeneracy of the energy lev- 
els. The dependence of the lifetimes of the resonances on 
the angular momentum and in particular the existence of 
so-called quasibound states has been analyzed in detail. 
Quantum states of neutral fermions and bosons in a mag- 
netic quadrupole configuration have been explored jl86l ] : 
The distribution of the energies and decay widths of the 
resonances were studied identifying the conditions under 
which states with long lifetimes can be achieved. A pe- 
culiar class of short-lived negative energy resonances in a 
magnetic quadrupole field has been found recently [1871 ] 
and a mapping of the two branches of positive and nega- 
tive energy resonances was derived. Global properties of 
the resonant quantum dynamics of neutral spin-1 atoms 
in a magnetic guide have been studied [188] including 
the effect of a Ioffc field. For certain parameter regimes 
the ground state resonance was shown to exhibit a longer 
lifetime than the energetically neighbored excited states. 
The case of a comparable h yper fine and field interaction 
was explored very recently [1891 ] . 

The past few years have seen substantial progress to- 
ward the development of stable trapping and guiding 
mechanisms and trapping on microscale levels, as in opti- 
cal lattices and on atom chips. There are two domains of 
interest: The case of a highly excited CM motion of the 
Rydberg atoms which can be described in terms of clas- 
sical dynamics or semiclassical approaches and the case 
of a quantized ground and excited CM motion. We will 
review in the remaining part of this section the status 
in case of a semi-classical CM motion (see in particular 
the review in Ref. 26]). The subsequent section will be 
devoted to the quantum Rydberg regime i.e. the case of 
a coupled quantum CM and electronic Rydberg motion. 

Trapping of long-lived strongly magnetized Rydberg 
atoms being collected in a superconducting magnetic trap 
with a strong bias field (2.9 T) has been explored ex- 
perimentally [3l|. Here so-called drift or guiding cen- 
ter Rydberg atoms (see subsection IIII F|) are produced 
by Rydberg collisions. Observation of oscillatory mo- 



tion of Rydberg atoms in the resulting magnetic poten- 
tial has provided strong evidence that long-lived Ryd- 
berg atoms have been trapped. Trapping perio ds u p to 
200 ms have been found. In a follow-up work [l68j the 
picture of the guiding center drift atoms in a strong mag- 
netic field, based on a hierarchy of interactions and time 
scales, has been used to derive simple models of the effec- 
tive magnetic moment and to p e rform classical trajectory 
calculations (see also refs. [l69l Il90j ] for the background 
on circular Rydberg atoms in strong magnetic fields). 
An anisotropic response of hydrogen atoms to inhomo- 
geneities longitudinal and transverse to the magnetic field 
has been observed. A proposition to use the pondcromo- 
tive energy of Rydberg electrons in standing-wave light 
fields to form an optic al lat tice for Rydberg atoms has 
been worked out in ref. 174]. Atom-chip technique pro- 
posals have recently been worked out for the trapping of 
a single Rydberg atom in a circular state [175| . The small 
size of microfabricated structures allows for trap geome- 
tries with microwave cut-off frequencies high enough to 
inhibit the spontaneous emission of the Rydberg atom, 
paving the way to complete control of both external and 
internal degrees of freedom over very long times. 



A. Atomic structure and quantum dynamics in an 
inhomogeneous magnetic field 

Our objective is to develop an ab initio approach of 
moving Rydberg atoms in inhomogeneous magnetic field 
configurations. We thereby provide a consistent descrip- 
tion of the quantum behaviour for the center of mass and 
electronic motions: both types of motion will be treated 
on equal footing. A major goal is to work out and an- 
alyze configurations which allow for the controlled pro- 
cessing of Rydberg atoms in the regime where both the 
CM and electronic motion are governed by quantum ef- 
fects. This includes traps and waveguides for ultracold 
Rydberg atoms. 

As we have seen in chapter IIII1 the case of a moving 
atom in a homogeneous magnetic field exhibits a number 
of intriguing phenomena. In particular, we encounter 
an intricate coupling of the CM and electronic motion 
which gives rise to both new structural properties and 
dynamics. It has to be expected that moving atoms in 
an inhomogeneous field will exhibit their own unique be- 
havior. In view of the high dimensionality and large den- 
sity of levels, it is advisable to fi r stly study the elect ronic 
structure for fixed nucleus [HO, EM, EM EM EH and 
subsequently investi gate the compl e te problem of a mov- 
ing Rydberg atom [196|, EM EM Ell HM HM US- 
Concerning the magnetic field config urat i ons, we will fo - 
cus here on the 3D quadrupole field [l~9ll [192 . fl9l fl95| . 
For the 2D quadrupole field (side guide) and the Ioffe- 
Pritchard configurations, we refer the reader to the liter- 
ature EH. 
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1. Electronic structure of Rydberg atoms for fixed nucleus 

Keeping the position of the nucleus fixed and study- 
ing exclusively the electronic structure of the atom in 
an inhomogeneous external magnetic field is certainly a 
special but instructive case of a Rydberg atom in a mag- 
netic field. It corresponds to the situation where the con- 
finement length of the CM motion is much shorter than 
the extension of the electronic Rydberg wave function. 
One of the motivations to study the impact of the inho- 
mogeneity of the field on the electronic structure is the 
experimental availability of extremely strong gradients; 
current carrying wire s on atom chips can reach gradients 
up to 10 4 T/cm [2031 ]. It is thus expected that such large 
field gradients have sizeable impact on the Rydberg wave 
functions. 

The 3D magnetic quadrupole field [IM EM EM Ell 
is characterized by B(r) = b(x, y, — 2z) with the gradient 
b. This vector field is rotationally symmetric around the 
z-axis and invariant under the z-parity operation. A cor- 
responding vector potential in the Coulomb gauge reads 
A(r) = | [B(r) x r]. Adopting atomic units (b = la.u, = 
4.44181 x 10 15 ^) the nonrelativistic spinor Hamiltonian 
in case of a single active Rydberg electron reads 



Hfn 



^ - 1 

2 y/ x 2 + y 2 + z 2 



b zL 2 



b 2 b 
+ —z 2 [x 2 + y 2 ) + - (a x x + o y y - 2<r z z) (56) 



xdy) and a x , cr y , a z are the Pauli 
\<r). The paramagnetic (oc b) or 



with L z = —i (yd. 
spin matrices (S 

orbital Zeeman-term depends, in contrast to the situa- 
tion of the atom in a homogeneous field, not only on L z 
but also linearly on the z-coordinate. The diamagnetic 
term (oc b 2 ) represents a quartic oscillator coupling term 
between the cylindrical coordinates p = \/ x 2 + y 2 and 
z. In a homogeneous field the diamagnetic interaction is 
a pure harmonic oscillator term proportional to p 2 and 
yields a confinement perpendicular to the magnetic field. 
The coupling of the electronic spin to the quadrupole field 
depends linearly on the Cartesian coordinates. The latter 
prevents the factorization of the motion in the spin and 
spatial degrees of freedom and renders the correspond- 
ing Schrodinger equation into a spinor equation. This 
is again in contrast to the case for a homogeneous mag- 
netic field where the spin component along the field is a 
conserved quantity. 

As a consequence of the rotational invariance of the 
quadrupole field the z-component of the total angular 
momentum J z = L z + S z is conserved, i.e. [Hf n , J z ] = 0. 
Additionally we have the discrete symmetry represented 
by the unitary operator P^OP z , i.e. [Hf n , P^OP z \ = 0. 
Here O = a x exchanges the components of a i-spinor and 
P<t, represents the 'reflection' P^ : <f> — > 2tt— <fi. Apart from 
these symmetries the Hamiltonian possesses two general- 
ized anti-unitary time reversal symmetries namely TOP z 
and TP§ involving the conventional time reversal opera- 
tor T (T 2 = 1). The operators TOP z , TP,p and P^OP z 



form an invariant Abelian sub-group. Together with J z 
they obey the following (anti-)commutation rules: 

[J Z ,TP^] = {J z ,TOP z } = {J M ,P^OP,} = (57) 
[TP^TOP,] = [TP^,P^OP z ] = [P4,OP z ,TOP z \ = (58) 

Apparently the spin-spatial symmetry operations form 
a non- Abelian symmetry group. This group generated by 
P<pOP z and J z is isomorphic to ® C s . 

The interplay of the above symmetries leads to a de- 
generacy of the eigenvalues of the Hamiltonian (f56|) . In- 
deed, if we consider a state \E, m) which is an eigenstate 
of the Hamiltonian (p)6"|) with the energy E and to J z 
with the half-integer quantum number m, then since the 
state TOP z \E,m) is also an energy eigenstate with the 
energy E, we obtain 

J z TOP z \E,m) = -TOP z J z \E,m) = -mTOP z \E,m) 

Thus the state TOP z \E, m) can be identified with 
\E, — m). Hence the states with the eigenvalues m and 
— m are degenerate. This two-fold degeneracy of each en- 
ergy level in the presence of the inhomogeneous magnetic 
field is a remarkable feature reminiscent of the Kramers 
degen eracy of spin | systems in the absence of external 
fields 



204] 



Let us now discuss on the basis of some examples new 
features that occur for the electronic structure of Ryd- 
berg atoms in inhomogeneous magnetic fields using the 
specific example of the 3d quadrupole field. To demon- 
strate the deformation effects on the electronic cloud we 
consider for reasons of illustration an extremely high gra- 
dient. Figure [T9l shows the spatial probability density of 
the ground state in the quadrupole field $Wn) at b = 10 
and the homogeneous field (fH)b ) at B = 10 (m = h 
in both cases). For the quadrupole field we observe an 
asymmetric deformation with respect to the 9 = ^-plane: 
the electronic wavefunction is almost completely confined 
to the upper half- volume (6 < ^) which is a consequence 
of the symmetries: z-parity is not conserved and eigen- 
states appear in pairs one being the mirror image of 
the other with respect to the reflection at the x — y- 
plane. Furthermore we observe that the electronic mo- 
tion is localized particularly along the two 'channels' for 
6 = corresponding to the lower z-axis and 9 = 5 being 
the x — y-plane. This property as well as the detailed 
shape of the electronic probability density in the individ- 
ual half-volume is determined by the diamagnetic term 
which is dominant in the high gradient regime. For the 
quadrupole field it is proportional to sin 2 6 cos 2 6 reach- 
ing its maximum value at 9 = ? j ^ • The probability 
density in a homogeneous field (see figure \T9b) exhibits 
the above-mentioned corresponding reflection symmetry 
due to the invariance of the corresponding Hamiltonian 
with respect to z-parity. Here we observe the maximum 
of the probability density at 9 = ? and a deformation 
towards 9 = ir and 9 = leading to a cigar-like shape. 
The diamagnetic term is proportional to sin 2 9 having 



its maximum value at 9 



thus coinciding with the 
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FIG. 19: Spatial probability density of the ground state in 
the quadrupole field (a) and a homogeneous field pointing in 
the ^-direction (b) for b = B = lOa.u. and m = |. From 

PH. 

regions possessing the strongest deformation of the prob- 
ability density. For both field configurations the proba- 
bility density vanishes at r = 0. 

To explore the compression of Rydberg states in strong 
(atom-chip) laboratory gradients we first remark that 

(r) = {r)TOP z = (r) Jz - (59) 

where the lower index at the expectation values indi- 
cates the type of eigenstates employed. For the one- 
electron problem without external field we have (r) H = 
h (in 2 — 1(1 + 1)). Since I satisfies < I < n — 1 one 
obtains corresponding upper and lower bounds for the 
expectation value of r: n (n + i) < (r) H < \n 2 . In fig- 
ure[2D]this expectation value is shown for the atom in the 
quadrupole field as a function of n for two different gra- 
dients (n serves here as a label for the energy levels and 
is not a good quantum number!). The radial expectation 
values lie in between the boundaries given by the field- 
free inequality indicated by the dashed lines (see figure 
^Uk ) . Expectation values of states belonging to the same 
n-multiplet are arranged in vertical lines expressing their 
energetical degeneracy (see figure l2"0"k for b — 10 -10 ). 
Here states with large expectation values of r possess 
small expectation values of the angular momentum and 




FIG. 20: Expectation value of the radial coordinate r plotted 
agains t n at different gradients (a: b — 10 -10 , b: b = 10 -8 ). 
From [HI]. 

vice versa. The situation is different for b — 10~ 8 . Here 
a systematic decrease of the r-expectation values takes 
place for states with an energy corresponding to n > 30. 
This is also the energy regime where the inter n-manifold 
mixing due to the diamagnetic term sets in. 

As a next step let us consider properties of the Ry- 
dberg atoms that involve the spin degrees of freedom. 
In a homogeneous magnetic field, the projection of the 
spin operator onto the direction of the magnetic field is 
a conserved quantity. In this case one can choose the 
energy eigenstates to be also eigenfunctions with respect 
to S z . In the quadrupole field, S z is not conserved, and 
we therefore consider the expectation value 

(S.)j, = ~ [(« I u) (d | d)] . (60) 

Figure |2"T1 shows the distribution of (S z ) j for electronic 
states of the m = h subspace as a function of n. Since 
S z is not conserved, the values of (S z ) are allowed to 
cover the complete interval [— s, 5]. For b = 10~ 10 (fig- 
ure ETJl) and a low degree of excitation the expectation 
values are evenly distributed over the interval. When 
reaching highly excited states this pattern becomes in- 
creasingly distorted. The expectation values agglomer- 
ate at —0.35 and 0.35 for n > 50. Due to the approxi- 
mate degeneracy of the energy levels at low energies, i.e. 
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FIG. 21: Expectation values of the z-component of the spin 
operator as a function of the quantum numbe r n f or different 
gradients (a: b = 10" 10 , b: b = 10" 8 ). From [l92T ]. 

small n, the values of (S^) form vertical lines. These lines 
widen for higher n. Since at b = 1CP 10 no significant n- 
mixing up to our maximum converged energy levels takes 
place, neighboring lines are well separated. For a higher 
gradient b — 1CP 8 (figure l2~Ib) the above properties are 
equally present for low-lying states. However, with in- 
creasing excitation energy we now observe a complete n- 
mixing regime where the regular patterns disappear and 
we obtain an irregular distribution of (S z ). Overall the 
distribution narrows, e.g. for n = 40 the occupied inter- 
val is approximately [—0.3,0.3]. We remark that due to 
the fact that the ^-operator anti-commutes with TOP Zl 
we have (S z )^ OP = 0. Apparently there is no preferred 
direction for the electronic spin in a state obeying the 
TOP z symmetry. 

Due to the coupling of the spatial and spin degrees of 
freedom the spin orientation becomes spatially dependent 
and it is instructive to consider the spin S^-polarization 
Ws- For a J z -eigenstate \E,m) it reads 

(£,m|r)S,(r|E,m) 
S[ > (E,m\r)(r |E,m) 

= l |Hr)M(rf|r)| 2 
2|( u |r)| 2 + |(d|r)f 

Figure [H] shows the spatial probability distribution 
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FIG. 22: Spatial probability density (a) and S'z-polarization 
(b) for the 1117th excited state for m = i and b = 1CF 8 . 
At large r the 5^ -polarization becomes similar to Wg indi- 
cating an antiparallel alignment of the electronic spin to the 
magnetic field. From [l92| ]. 



(T2"2"h ) and the ^-polarization (f2"2"b ) for the 1117th ex- 
cited state for m = | and b = 10 -8 In contrast to the 
constant ^-polarization we would encounter in the ab- 
sence of a field or a homogeneous field we observe a com- 
plex pattern of domains exhibiting different spin orien- 
tation (white: spin up, black: spin down). At low r 
values these domains form a pattern similar to that of 
a chess board. The junctions where four spin domains 
meet each other coincide with the nodes of the spatial 
probability density. The Coulomb interaction as well as 
the spin-Zeeman term are responsible for the interwoven 
network of island of different spin orientation. The ad- 
ditional presence of the orbital Zeeman and the diamag- 
netic term leads to a deformation of this network. For 
low radii the spin orientation changes locally from island 
to island generating an appealing pattern whereas we ob- 
serve an overall tendency of the electronic spin polariza- 
tion in the region characterized by stripes (r ~ 2000). 
Here, independently of the nodal structure, the spin ori- 
entation changes smoothly from downwards at 8 = n to 
upwards at 8 = 0. This feature can be understood by 
inspecting the spin Zeeman term only. 

In a homogeneous magnetic field parity is a symme- 
try and therefore atomic electronic eigenstates (for fixed 
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nucleus) do not possess a permanent electric dipole mo- 
ment. Let us investigate the electric dipole moment of 
the electronic states in the quadrupole field. First of all 
we remark that for dipole transitions the correspond- 
ing matrix element (E' , m! \ r sin e ^ \E, m) is only non- 
zero if m! — m = ±1. The expectation value of D a ± in 
the J z -eigenstates vanishes 

(D a ±) = (E,m\r sin (9 e ±i4 '\E,m) = (62) 

However, the expectation value of is in general non- 
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FIG. 23: Expectation value of the dipole operator D n plotted 
against n for different gradients (a: b = 10 -10 , b: b = 10 -8 ). 
From fl92| . 

zero. (D n ) is shown in figure I23f a.bl for the two gradi- 
ents b = 1CT 10 and b = 1(T 8 , respectively. For b = lCT 10 
the electric dipole moments belonging to the same n- 
multiplet are arranged along vertical lines, which is a re- 
sult of the approximate degeneracy of the energy levels. 
The substates for fixed m of a given n-multiplet exhibit 
different dipole momenta spreading between two (upper 
and lower) bounds that depend linearly on n. With in- 
creasing gradient and degree of excitation the n-mixing 
starts and disturbs the observed regular pattern. For 
b = 10~ 8 and n > 35, the distribution of the dipole mo- 
ments becomes completely irregular. Therefore we en- 
counter the remarkable effect that the external magnetic 
quadrupole field induces a state dependent permanent 



electric dipole moment. This is the result of the asym- 
metric form of the wavefunction in the quadrupole field 
(see discussion of the probability density in a strong gra- 
dient field above). 



2. Rydberg atoms motion 

a. 3D quadrupole trap: We will start with some con- 
ceptual ideas and remarks on the 3D q uadrupole case , 
refering further details to the literature jl96l . Il97l Il98j , 
and will then work out the Ioffc-Pritchard configuration 
in more detail. Finally the addition of an electric field, 
i.e. a magneto-electric trap, is outlined. 

Similar to the case for a fixed nucleus (see subsection 
IIV A 1|) , a moving Rydberg atom in a 3D quadrupole field 
possesses a rotational symmetry around the z-axis. The 
corresponding conserved quantity is the projection J z of 
the total (spin and orbital) angular momenta onto the 
symmetry axis of the quadrupole field. To exploit this 
constant of motion the re a re essentially two pathways 
that could be pursued [l97l | in order to derive effective 
equations of motion. In both cases the first step is to 
transform the two-body problem from laboratory to CM 
and internal electronic coordinates and to simplify the 
resulting Hamiltonian and its coupling terms with a suit- 
able Unitarian. The first approach eliminates two of the 
six (CM plus electronic) spatial degrees of freedom and 
leads to an (infinite) set of coupled channel equations for 
the spin and spatial degrees of freedom. The second ap- 
proach introduces J z as a canonical momentum, thereby 
eliminating the corresponding cyclic a ngle. Th e latter 
approach has subsequently been used jl96l . Il98j | to ex- 
plore the case of a 3D quadrupole field. Since many of 
the concepts and ideas developed will be described below 
for the case of the Ioffe-Pritchard configuration we pro- 
vide here only the key result for the 3D quadrupole field. 
It turns out that trapping is possible for a certain class 
of coupled CM and electronic quantum states: Both the 
CM and electronic motions have equally large angular 
momentum states. The requirement that specifically the 
CM motion bears a large angular momentum is due to the 
fact that the electronic interaction terms in the Hamil- 
tonian lead to a clustering of adiabatic potential energy 
surfaces and therefore strong nonadiabatic effects. The 
latter, however, leads to an uncontrollable unstable CM 
motion. Consequently a spatial regime where the energy 
is dominated by the electronic energy has to be avoided 
altogether and this is achieved by increasing the angular 
momentum which leads to an increase of the angular mo- 
mentum barrier terms which in turn shrinks the spatial 
regime where the electronic interaction exerts influence 
on the CM motion. 

To summarize, trapping of Rydberg atoms in the 3D 
quadrupole field is possible, but only for a class of high- 
angular momentum CM and electronic states, while con- 
trol and manipulation of these states are difficult. Of 
course, the experimental preparation of high angular mo- 
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mentum CM states requires additional effort. As indi- 
cated the electronic part of the trapped states possesses 
also a large angular momentum i.e. they are close to 
circular states [l| for which well-defined experimental 
recipes of preparation exist [205l ] . 

b. Ioffe-Pritchard trap: The starting-point is the 
two-body non-relativistic Hamiltonian in the laboratory 
frame, which incorporates the Rydberg electron and the 
nucle us or effec tive (closed shell) ionic core of an alkali 
atom (l99l [202l. 

Hintt = ^(Pi-9i A ( r i)) 2 + ^(P2-q 2 A(r 2 )) 2 

+V{\r x - r 2 |) - mi • B(n) - fi 2 ■ B(r 2 ) .(63) 

where fii, H2 are the magnetic moments of the particles 
1 and 2, respectively. First one introduces CM and elec- 
tronic relative coordinates and performs a unitary trans- 
formation reminescent of the transformation applied in 
the case of a homogeneous magnetic field (see U n in Eq. 
(|18p ). Retaining only leading order terms with respect 
to the inverse masses, neglecting the diamagnetic interac- 
tion, and inserting the Ioffe-Pritchard field configuration 
(up to the linear order) for the vector potential, respec- 
tively, give 




(64) 
(65) 



while the Hamiltonian in a.u. is 

P 2 1 

2M + 2 



H 



ip 



BL Z + G(x + X)(y + Y)p z 



-mB(R + r)-ii 7 B{R) (66) 



where Ha is the operator for the field-free atom. This 
allows for the separation of the free CM motion in Z- 
direction. The Zeeman term comes from the uniform 
Ioffe field generated by the Helmholtz coils. The next 
term, involving the field gradient G, arises from the lin- 
ear field generated by the Ioffe bars and couples the rel- 
ative and CM dynamics. The last term couples the spin 
of particle 2 to the magnetic field. Since the electronic 
spins of closed shells combine to zero, the spin of parti- 
cle two is the nuclear spin only. Even though fi 2 B scales 
with I/M2, this term is important for a proper symmetry 
analysis, but will be omitted later on. 

The above Hamiltonian Hjp is invariant under a num- 
ber of symmetry transformations Us that are composed 
of the elementary operations listed in Tab.|TJ The unitary 
symmetries are 



PxPySz^z 

ByPz^xy^xy^xy 

PPT Q* V* 
r x r z'-xy&xy'-'xy • 



(67a) 
(67b) 
(67c) 



TABLE I: Symmetry operation nomenclature. Pj, Sj, and Ej 
x, but hold of course also for j = y, z. 



are exemplified by j 



operator 




operation 




P x 


x parity 


x — > —x, X —f 


-X 


Sx 


electronic spin x op. 


Sy ► Sy , Sz 


-+ -Sx 




nuclear spin x op. 


E H — > —Ey, E z 


-E, 


Ixy 


coordinate exchange 


x <-> y, X <-> Y 




Sxy 


el. spin component exc. 


S x * Sy , Sy 




^>xy 


nuclear spin comp. exc. 


Ex ^ ^-*y 1 




T 


conventional time reversal 







The Hamiltonian is also left invariant under the antiu- 
nitary symmetry transformation TP y . By consecutively 
applying the latter operator and the unitary operators 
Eqs. (|67a|) . (|67b|) and (|67c[) . it is possible to create fur- 
ther antiunitary symmetries: 



TP x Sz^z 

T Pz^xySxy^xy 

T 1 p PPT Q* V* 
J r x r y r z 1 xy°xy L -'xy 



(68a) 
(68b) 
(68c) 



Considering that S% y = — S z and Y? xy = — E 2 and that 
T neither commutes with S y nor with S xy and T, xy , 
one finds that the operators (|67aH68c[) form a symmetry 
group. 

If no Ioffe field is present (P> = 0), eight additional 
symmetries can be found leaving the Hamiltonian invari- 
ant. For an effective one particle approach (and the cor- 
responding one particle symmetries) this was discussed 
in Ref. fl93 |. 

For the parameter regime we are focusing on (typical 
experimental Ioffe field strength and gradients) there is 
no inter-manifold mixing and each n-manifold can there- 
fore be separately considered. We introduce the scaled 
CM coordinates, R — > 7-3R, with 7 = GM, and the 
scaled energy e = 73/M to remove explicitly the de- 
pendence on B, G, and M, Introducing the effective 
magnetic field 



G(X,Y) 




C = BM-y- 



(69) 



and omitting the zero-field energy offset E\, the Hamil- 
tonian matrix within an 71-manifold has the simple form 



n 



pi 



p2 

y 



+li-G(X,Y)+^(xy Pz +xS x -yS y ). (70) 



The first term is the CM kinetic energy, fi is the 2n 2 - 
dimensional matrix representation of the total magnetic 
moment of the electron, |(i r + 2S), and the second term 
in eq. (|70"|) describes its coupling to the effective magnetic 
field G. 

The large differences in masses and velocities in the 
two body system allows for an adiabatic separation of the 
electronic and the CM motions. This is true regardless of 
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the slow motion of the Rydberg electron compared to the 
ground state electron. It remains to be investigated how 
the spatial field inhomogeneity will affect density of lev- 
els and level crossings. The goal is therefore to find con- 
ditions under which isolated adiabatic potential energy 
surfaces might exist. The adiabaticity will be provided 
by the CM motion. 

The adiabatic approximation is introduced by sub- 
tracting the transversal CM kinetic energy, T — 
{Pi + Py)/2, from the total Hamiltonian fl70j!. The re- 
maining electronic Hamiltonian for fixed position of the 
CM reads 

He = ti-G{X,Y)+^{xyp z +xS x -yS y ). (71) 

The electronic wave function ip. K depends parametrically 
on R and the total atomic wavefunction can be written 
as 

Mr,R)) =\<p K (r;R))®\MR)) , (72) 

where \tj) v (R)) is the CM wave function, and the elec- 
tronic wave functions are solutions of 

H e \<p K (r; R)) = E K (X, Y) |^(r; R)) (73) 

akin to the electronic wave functions in the Born- 
Oppenheimer approximation. We emphasize that this 
adiabatic separation of the CM and electronic coordi- 
nates for the Rdyberg atom in intrinsically tied to the 
field inhomogeneity. The electronic potential energy sur- 
faces on which the CM moves are E K (X, Y). The equa- 
tion of motion for the CM wave function now reads 

(T + E K (X, Y)) \ip„{R)) = e v \MR)) ■ (74) 

A discussion of the nonadiabatic coupling terms can be 
found in Ref. (202^ . 

It can be demonstrated that the energy surfaces 
E K , exhibit three mirror symmetries. If we apply 
the corresponding restricted symmetry operation Up = 
P x P y S z T, z in Eq. (|67a[) . that was already shown to leave 
the Ioffc-Pritchard Hamiltonian in Eq. (f66|) invariant, to 
the (unsealed) electronic Hamiltonian H e , we find 

tf p H e (r;X,Y)Up = H e (r;-X,-Y) . (75) 

i.e. the energy surfaces are inversion symmetric. The 
symmetry operator Uy — TP V , and the operator that is 
composed of Uy and Up, namely Ux = TP x S z ^ Zl mirror 
the energy surfaces at the axes, 

U Y H e (r;X,Y)U Y = H e (r;X,-Y) , (76) 
U x H e (r;X,Y)U x =H e (r;-X,Y) . (77) 

The electronic Hamiltonian in Eq. (|T3[> . with the core 
fixed at an arbitrary position, is three-dimensional. No 
symmetry arguments can be exploited to reduce its di- 
mensionality. 



Let us now analyze the properties of the electronic adi- 
abatic potential energy surfaces resulting from the di- 
agonalization of the electronic Hamiltonian for different 
regimes of the Ioffc field strengths and field gradients. 
These two parameters can be employed to change the 
appearance of the potential in which the center of mass 
dynamics takes place. A careful inspection of the elec- 
tronic Hamiltonian shows that its first term, fi ■ G(X, Y), 
includes the terms X(^L X + S x ) — Y(^L y + S y ), that are 
of the order of (Li) rs n for high angular momentum 
states, and the Zeeman term C,(\L Z + S z ), which can be 
as large as (n. The second term, "f x ^(xyp z + xS x — yS y ) 
can be estimated as 7 1 / 3 n 3 . In a nutshell, we have the 
following relative orders of magnitude, 

1 , C and 7^ n 2 . (78) 

Due to the special form of the electronic Hamiltonian, 
changing the magnetic field parameters B and G while 
keeping their ratio C/l 1 ^ — B/G (and n) constant re- 
sults in a mere scaling of the CM coordinates (we remark 
that the characteristic length scale of the CM dynamics 
is of the order of one in scaled atomic units). 

To understand the impact of the Ioffe field strength B 
on the adiabatic energy surfaces, we isolate its effect by 
suppressing other influences. This can be done by choos- 
ing a relatively low field gradient G and/or a compara- 
tively low n. such that ■j'Sn 2 is small, and the last term 
in Eq. (|71[) will become negligible. Within this regime 
approximate analytical expressions for the electronic adi- 
abatic energy surfaces can be derived. We therefore di- 
agonalize the resulting electronic Hamiltonian 

H e = \G(L + 25) . (79) 

by applying the spatially dependent unitary transforma- 
tion 

U D (X,Y) = e i<l>(L*+S») e i0(L y +s y ) ) ( 8Q ) 

with (j> = arctan^, cos/3 = M 2 B\G(X, Y)\~ x and 
sin/3 = -\/A 2 + Y 2 \G(X,Y)\-\ This yields 

ulH e U D = ^(L z + 2S z )\G(X,Y)\ (81) 

for the transformed approximate electronic Hamiltonian. 
The spatially dependent transformation U d locally ro- 
tates the magnetic moment of the electron, which in- 
cludes its spin and its angular momentum, such that it is 
parallel to the local direction of the magnetic field. The 
operators L z and S z are not identical to the ones before 
the transformation in Eq. (|80| ; they are rather related to 
the local quantization axis defined by the local magnetic 
field direction [l96j . 

The approximate adiabatic potential surfaces are now 

E K (X,Y) = i(m, + 2m s )\G(X,Y)\ 

= i(m* + 2m s )vOs: 2 +y 2 +C 2 ■ (82) 
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The possible combinations of mi and m s yield 2n + 1 
energy surfaces. The surfaces highest and lowest in en- 
ergy correspond to circular states, (\mi\ = l m ax = n — 1, 
mi + 2m s = in), and they are the only non-degenerate 
ones. For the other surfaces (|m/ + 2m s | < n), by 
2n — \mi + 2m s + 1| — \mi + 2rn s — 1|. Starting from 
the highest energy surface, the levels of degeneracy are 
thus 0, 2, 4, 6, ... . 

The approximate surfaces E K are rotationally symmet- 
ric around the z-axis. An expansion around this axis 
(p = v 'X 2 + Y 2 -c C) yields a harmonic potential, 



-3 3 -3 3 



1 



1 



E K (p)^((+—p 2 )--(mi + 2m s ), 
while we find a linear behavior, 



(83) 



E K (p) ^ ^ ■ (mi + 2m s ) , (84) 
when the center of mass is far from the z-axis (p 3> £)■ 
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FIG. 24: Sections along the X-axis through the electronic 
adiabatic energy surfaces of an entire n = 3 manifold. The 
field gradient is fixed at G — 1 Tesla/m in order to suppress 
the influence of the last term in H e (Eq. I71[) . From left to 
right, t_= BMy~ 2 ^ 3 increases due to an increasing Iofle field. 
From [203. 



Fig. [M] shows sections (7 



l/3 n 2 = 



0.003, c = 0.01 - 1) 



through all the surfaces n = 3. This principal quantum 
number has been chosen in order to keep the sections 
simple while displaying the entire n-manifold. 

Fig. [25] shows for n = 30, G = 0.1 T/m (-> 7 1 / 3 ri 2 = 
0.14) the uppermost 21 energy surfaces for different val- 
ues of B. The harmonic behavior around the origin 
is clearly visible and the minimal distance between the 
surfaces becomes larger with increasing £. Since C, and 



7 



l/3_2 



are of the same order of magnitude in subfigure 
(a), the contribution of the last term in Eq. (j7T|) . that 
lifts the degeneracy of the curves, is visible. 

The minimum energy gap between two adjacent sur- 
faces is at origin 
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FIG. 25: Sections along the X-axis through the uppermost 
21 surfaces of the n = 30 manifold of 87 Rb for increasing 
ratios B/(Gn 2 ). The field gradient is fixed at G = 10 T/m 
while the Ioffe field is increased from top left to bottom right. 
(B = 24 mG, B = 48 mG, B = 0.24 G, B = 0.48 G). For 
small ratios B /(Gn 2 ), the influence of the second term in Eq. 
(|71[) is not completely suppressed as can be see n fro m the 
lifted degeneracies in the upper subfigures. From [202f | . 
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FIG. 26: Section through the n = 30 manifold for a field 
strength of 0.01 Gauss and a field gradient of 20 T/m ( 87 Rb). 
A large number of avoided crossings can be observed. The up- 
permost curve, however, stays isolated from the other curves. 
The insets show the li near behavior of the surfaces far away 
from the .z-axis. From [2021 • 



(85) 



The parameter £ (and hence the field strength B) is 
the tool to control the energetic distance between the 
adiabatic surfaces. The uppermost energy surface proves 
to be very robust: It is energetically well-isolated from 
the other adiabatic surfaces and is therefore best-suited 
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for the trapping of the ultracold Rydberg atoms. This 
holds even in the regime of small Ioffe field strength B 
and comparatively large gradients G, illustrated by the 
corresponding sections of the adiabatic electronic energy 
surfaces shown in Fig. [26] for a Ioffe field strength of 
0.01 G and a field gradient of 20 T/m. For these pa- 
rameters, the contributions of all terms in the electronic 
Hamiltonian are of the same order of magnitude around 
X — 1. One immediate observation is that a large num- 
ber of avoided crossings occur between the surfaces, while 
the uppermost curve however remains isolated from the 
rest of the curves. Far away from the trap center, i.e. for 
large p = \/X 2 + Y 2 , the coupling term in Eq. (fTTj) . 
X{\L X + S x ) — Y{\L y ± S y ), becomes dominant. 

The energetically uppermost adiabatic electronic en- 
ergy surface is the most suitable for achieving confine- 
ment. It does not suffer significant deformation when 
the field gradient is increased and it stays energetically 
isolated from lower surfaces for a wide range of parame- 
ters. The motion along the Z direction is unrestricted. 

By choosing large gradients and appropriate bias fields, 
tight confinement for highly excited atoms can be ob- 
tained. For a Ioffe field strength of B = 0.1 G and a 
field gradient of G = 100 T/m, for instance, the ratio of 
(p) and (r) for the ground state (u = 1) is as small as 
(p) I (r) = 0.4. The extension of the CM wave function is 
thus smaller than the extension of the electronic cloud, 
and the Rydberg atoms can no longer be considered as 
point-like particles. Moreover, it is possible to enter a 
regime where the C M an d the electronic wave functions 
do not even overlap [202j |. 

We conclude that the Ioffe-Pritchard field configura- 
tion provides a strong confinement for ultracold Rydberg 
atoms in two dimensions. A relatively weak longitudinal 
confinement along the z-axis could additionally be pro- 
vided by employing a non-Hclmholtz configuration. The 
preparation of a quasi-stable ID Rydberg gas in the wave 
guide provided by the Ioffe-Pritchard field configuration 
has been worked out very recently \20l\ by adding a ho- 
mogeneous electric field component perpendicular to the 
Ioffe field. The undesirable collisional autoionization pro- 
cess of ultracold Rydberg atoms is prevented here by the 
repulsive dipolc-dipole interaction among the atoms. 



B. Lifetime of trapped Rydberg atoms 

Thus far, we were mainly concerned with static proper- 
ties of Rydberg atoms in strong inhomogeneous magnetic 
fields. However, when confined in magnetic traps for suf- 
ficiently long times, spontaneous decay of excited levels 
starts to play a crucial role for their dynamics. Owing 
to the large excursions of highly excited electrons, the 
Rydberg state decay rate can be small, scaling in the 
magnetic field-free case as [l[ 

r n i m ~ n~ 3 (86) 
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FIG. 27: Spontaneous decay rate of circular Rydberg states 
with m — ±39 as function of the magnetic field. Exact numer- 
ical results (dots) are compared to the prediction of Eqs. (187[) 
and (|88[) (solid lines). The upper x-axes shows the relative 
Zeeman shift Ae = \m\pBB / E n i m . 

with n for a fixed I. The decay rate of circular states with 
I ~ n — 1 has an even stronger n-dependence, ~ n~ 5 , re- 
sulting in extremely long lifetimes of several milliseconds 
for typical excitations of n > 40. 

Ext ernal magnetic fields strongly affect spontaneous 
decay fl69l 12061. [207j| . and, depending on to, can signifi- 
cantly enhance or reduce the lifetime of Rydberg atoms. 
To illustrate this effect, we consider first the simple limit 
of weak B, in which the binding energies experience a lin- 
ear Zeeman-shift E n i m (B) — E n i m (0) + mp^B, while the 
transition dipole matrix elements remain unaltered. Here 
Pb = eh/2m denotes the Bohr magneton. Since the spon- 
taneous decay rate is proportional to the third power of 
the energy difference between the initial and final states, 
we obtain for the decay from (n, l,m) to (n',l',m ± 1), 
the magnetic field dependence 

r nlm (r>\ — fl -r _"^"^MB_B\ r nlm fn\ 
1 n'l'm±l \ B ) ~ y 1 T n2 _ n , 2 J I n'l' m ±l W , 

(87) 

where 1Z = 13.6 eV. For circular Rydberg states (|m| = 
I = n— 1), the transition (n,m,l) — > (n — 1, m± 1, 1) pro- 
vides the only dipole-allowed decay channel. The decay 
rate r* rc of circular R ydbe rg states with m — ±|to|, as 
obtained from eq.([57|) |206j l 

^)=(l±^f ^) 3 IW0) (88) 

shows that the decay of positive-m states is enhanced and 
the decay of negative-m states is reduced by an external 
magnetic field as compared to the field-free decay rate 
r circ (0). In Fig. [23 we compare these simple expres- 
sions to exactly calculated decay rates for circular Ryd- 
berg states with n = 40 and m = ±39. Eqs. (|87p and 
([55)1 yield a rather good description for magnetic fields 
smaller than ~ IT, for which the unperturbed energies 
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FIG. 28: Classical electron trajectories of strongly magnetized 
Rydberg atoms with negative (left) and positive (right) pro- 
jection of the angular momentum. In both cases, the trans- 
verse size is 10 4 a.u. The negative-m state shown in the left 
panel corresponds to a circular guiding center atom as dis- 
cussed in section IIII Fl From |l69l ] . 



are Zeeman shifted by up to 25% and the positive-m de- 
cay rates are enhanced by more than 50%. At larger 
fields, B > IT, Eqs. (J57J) and (JHHJ) systematically under- 
estimate the exact decay rates. This partially originates 
from the quadratic Zeeman shift, which increases the en- 
ergy difference between initial and final states for both 
negative and positive m states. Consequently, at fields 
that mark the strongly magnetized regime the above per- 
turbative treatment breaks down. 

The qualitative decay dynamics of such strongly mag- 
netized atoms can be understood within a classical pic- 
ture. Rydberg states with large negative L z fall into the 
class of circular guiding center atoms (see section IIII Fl) . 
Fig. |2"81 shows a typical trajectory of such an atom jl69l |. 
As will be discussed below, the atom's lifetime is limited 
by the decay of the E x B drift motion, which (see section 
IIII F[) has a frequency cjd = e/Bp 3 . The electron orbit 
produces an electric dipole that oscillates with frequency 
cjd and amplitude p. The drift motion, thus, radiates 
with the rate ~ u^p 2 = e 3 /B 3 p 7 . For comparison, a 
field-free, circular atom of the same size has a classical 
electron orbit with frequency ujq 



\J e 2 /mp 3 . Note that 



this expression coincides with the axial bounce frequency 
lj z , which we found to be much smaller than lou (see Eq. 
(|47|) . section HlIF[) . This simple estimate, hence, yields 
decay a rate r~ irc , suppressed by a factor of 



r cir c(o) 



mp 



,3/2 



B'> 



(89) 



The corresponding classical orbit with the same radius 
p but posit ive angular momentum L z is also shown in 
FigEE jl69l | . Flipping the angular momentum completely 
changes the character of the classical electron motion. 
The slow magnetron motion is replaced by a " giant" cy- 
lotron oscillation around the ion. Consequently, the de- 
cay rate r^ rc (i?) ~ ui 3 p 2 of circular, positive-m states is 
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FIG. 29: Left: m-dependence of the cyclotron, z-bounce and 
magnetron decay rates of Rydberg atoms in a 6T magnetic 
field. The rates are caluclated for the lowest states of the 
m manifolds for which the respective transition can occur. 
Right: Decay rate for the lowest states within the m- manifold. 
The adiabatic approximation (squares) is compared to ex- 
act results, showing agreement for |m| > 4£? -1 / 3 ~ 140 (see 
Eq. (|91|l '). The field- free decay rates of circular Rydberg atoms 
with a similar size demonst rates the strong suppression due 
to the magnetic field. From [l69l ]. 



drastically enhanced by a factor of 



3„9 



= (0) 



B 6 p 



(90) 



with respect to the field- free decay rate r c ; rc (0). 

An insightful quantum treatment of strongly magne- 
tized Rydberg atoms in the spirit of the g uidin g cen- 
ter approximation has been put forward in [l9dl ]. Us- 
ing a Born-Oppenheimer approximation, the z-bounce 
and t rans verse p dynamics are decoupled adiabatically 
[l69l [l90]. If this Born- Oppe nheimer approximation 



holds, which is the case for [169l | 
\m\ > AB~ 1/3 



(91) 



the energies associated with the cyclotron, the magnetron 
drift, and the axial bounce motion are adiabatic con- 
stants of motion, such that respective approximate quan- 
tum numbers can be defined. This permits to catego- 
rize the radiative decay in cyclotron transitions, bounce, 
and magnetron decays. Fig. [25] shows the correspond- 
ing decay rates as a function of m < 0. If the condi- 
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tion (f9Tj) is fulfilled, the three decay rates form a hier- 
archy, just like the corresponding oscillation frequencies 
(see Eq. (j47j) . section lHIFj) . Consequently a guiding cen- 
ter atom quickly dissipates its cyclotron energy, followed 
by a decay of the axial degree of freedom to its ground 
state, such that the overall decay rate of the Rydberg 
state is limited by the magnetron decay. As shown in 
Fig. [2H the latter can be suppressed by several orders of 
magnitude for highly excited states. 

Such guiding center atoms (with m < 0), when suf- 
ficiently deeply bound, typically have high-field seek- 
ing character, and, thus, are removed from the mag- 
netic traps. A trapped Rydberg gas, hence, consists 
of atoms in positive-m states, such that its sponta- 
neous decay is enhanced compared to the field-free case. 
However, the magnetic moment of strongly magnetized, 
negative-m states reverses its sign with increasing atom 
size (p > 2E cyc /e 2 ), such that very highly excited states 
can be trapped despite having m < 0. It is this sign 
reversal at high excitations that has enabled trapping of 
extremely long-lived, strongly-magnetized Rydberg gases 
[26L l3lj |. While this property is certainly appealing for 
studies of magnetized Rydberg gases, the opposite is the 
case for antihydrogen experiments, which aim at the pro- 
duction of ground state atoms through recombination in 
cold plasmas. It appears that the production of positive 
m is essential, since these states remain low-field seekers 
during their radiative cascade to lower-lying states. Here, 
the field-enhanced decay rates of these low-field seeking 
states is beneficial, as it accelerates the production of 
cold ground state antihydrogen atoms from radiatively 
decaying Rydberg gases. A more detailed account of the 
decay dynamics of such trapped Rydberg gases will be 
given in the next section. 



C. Long-time dynamics and cooling in strong field 
traps 

The preceding discussion of the magnetic moments and 
spontaneous decay rates fsection llVB[) elucidates how to 
separately treat either the translational center-of-mass 
dynamics in weak-gradient traps or the radiative cascade 
of Rydberg states. However, since the magnetic trapping 
potential depends on the atom's internal state, both of 
these processes are intima tely connected. In fact, it was 
shown in Refs. [206l |208| | that this entangled dynamics 
can lead to significant cooling of the trapped Rydberg 
gas during the radiative cascade of the atoms to their 
internal ground states. This effect may be of potential 
importance for antihydrogen experiments, as it strongly 
enhances the trapping efficiency of ground state atoms 
produced through radiative decay of initially recombined 
(see section [Vj highly excited Rydberg atoms. 

Several different magnetic and electric field config- 
urations are currently used to simultaneously confine 
charged and neutral particles, in order to synthesize neu- 
tral antihydrogen atoms within a two-component plasma 



and trap the produced atoms. Both the ATRAP and the 
ALPHA collaborations use a nested Penning trap (see 
section IPTCl) to create antihydrogen atoms. This plasma 
trap is su perimposed by either a quadrupole (ATRAP, 
[loa Hi^l) or an octupole (ALPHA, [67]) magnetic field 
to provide confinement for eventually formed neutral 
atoms. This additional field breaks the cylindrical sym- 
metry of the total magnetic field. A lack of cylindrical 
symmetry implies a lack of a confinement theorem (43j 
for the charged part icles and, thus, ultimately causes loss 
of plasma charges [21l| . However, both trapping con- 
figurations were experimentally demonst rated to p rovide 
sufficiently long confinement times [13, 12091 l21dj ]. Al- 
ternatively, the ASACUSA collaboration developed a so- 
called MCEO trap, which consists of an electric octupole 
field and a magnetic quadrupole (cusp) field to simul- 
taneously confine positrons, antiprot ons a nd antihydro- 
gen atoms in the same spatial region [212J. Producing a 
cylindrically symmetric magnetic field, this configuration 
promises long plasma confinement, but might cause atom 
loss at the magnetic field minimum where B = [2 12f | . 
In the following, we provide a general description of the 
dynamics of Rydberg gases in weak-gradient magnetic 
multipole traps. We will further assume that low-field 
seeking atoms, i.e. positive-m states, are present in the 
trap. The preceding formation of such Rydberg atoms in 
cold plasmas along with the characteristic internal states 
resulting from their formation will be discussed in section 

El . 

Since we consider slowly moving atoms in weak- 
gradient magnetic traps, it suffices to determine the 
atoms' field- and state-dependent magnetic moments 
n(B : r}) and radiative decay rates T 7 ^, (B). One can then 

employ a semiclassical treatmenfaiLthe gas dynamics de- 
scribed by the kinetic equation |206| 

77' 

-r^/(r,v,r,,t)(92) 

determining the time-evolution of the phase-space den- 
sity /(r,v,r/,t), i.e. the probability to find an atom 
at a position r, with a velocity v and in an atomic 
state 77, denoting a set of quantum numbers that defines 
the internal state. The trapping potential is given by 
U v (r) = fi(B,r])B(r). It is the ^-dependence of the con- 
finement potential that couples the translational center- 
of-mass dynamics to the internal state dynamics, as de- 
scribed by the right-hand side of Eq. ([92]) . It can be 
efficiently solved using a test-particle treatment. Within 
this approach, one randomly generates an ensemble of 
atomic states (r, v and 77) and propagates the classical 
coordinates according to the equations of motion. After 
each propagation time step, the internal state is changed 
according to the radiative decay rates by using a Monte 
Carlo procedure. 

Fig. [30] shows an example of the calculated time evo- 
lution of the atomic temperature and the fraction of 
trapped atoms in a magnetic cusp trap with a depth of 
AT. The temperature drastically decreases from its initial 
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FIG. 30: Time evolution of the number of trapped atoms 
(a) and the atomic temperature (b) for no = 44, mo = 43. 
Calculations have been performed for a magnetic cusp trap 
with a B-field difference of 4T. 
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value of 15K to below ~ IK, and the fraction of finally 
trapped ground state atoms is as large as 10%. Con- 
sidering that the effective trap depth decreases to about 
1/10 of the initial temperature, this is a remarkably large 
trapping efficiency, achieved only through the additional 
cooling. 

The basic mechanism can be understood in rather sim- 
ple terms. Whenever an atom at position ro decays to 
a state with a lower magnetic moment p! < /x, the local 
trapping potential decreases by a factor \J j [i. Conse- 
quently, the atom looses potential energy and, therefore, 
decreases its total energy by (//— /x)-B(r ). Subsequently, 
the total energy loss is distributed between potential and 
kinetic energy as the atom moves through the trap. This 
causes a net cooling effect and thereby increases the trap- 
ping efficiency of atoms following the radiative cascade. 

The cooling efficiency is generally expected to depend 
on the initial distribution of bound Rydberg states. How- 
ever, since these are generally unknown in experiments, 
calculations based on particular initial state distributions 
might be of limited practical applicability. Instead, Fig. 
13 II shows t he tr apping efficiency for particular initial en- 
ergies (no) 234] and angular momenta (too)- For fixed no 
and varying mo atom loss accelerates with decreasing mo , 
while the final fraction of trapped atoms is practically 
independent of mo. When plotted against the average 
principal quantum number (n) , the results collapse on a 
single curve. This insensitivity to Too is due to the fact, 
that the excited H is quickly driven to a c ircular Rydberg 
state during the radiative cascade [213j |. The following 
discussion can, thus, safely be restricted to circular Ry- 
dberg states. On the other hand, the initial value of n 
has a clear impact on the trapping efficiency, as shown 
in Fig |3"TT c) and l3"lTd) . The fraction of trapped atoms 
increases linearly with no, largely a consequence of dia- 
magnetic field effects on the atomic state. 

Further insights into the cooling process can be gained 
by considering a simplified spherically symmetric field 
B — /3(r/A) 7 . Here the parameter A controls the gradi- 
ent length and 7 the multipole order of the field, with 
7 = 1 corresponding to the cusp trap case, and larger 
values of 7 reproducing the central trap shape of higher 



FIG. 31: Fraction of trapped atoms as functions of time [(a), 
(c)] and average principal quantum number [(b), (d)], for no = 
44 and different initial 33 < m < 43 [in (a) and (b) at T(0) = 
10 K] and circular Rydberg atoms with 36 < no < 44 [in (c) 
and (d) at T(0) = 12 K]. Different curves are steps of two 
in mo and no. The inset shows the linear no-dependence of 
the fraction of trapped atoms (open circles), which disappears 
when strong magnetic field effects are neglected (filled circles) . 
The trap parameters are identical to those in Fig. 1301 From 



order multipole fields. Furthermore, we neglect the dia- 
magnetic term in the Hamiltonian, such that the exter- 
nal potential of circular states is given by U n — fi^nB. 
An atom in state n with a total center-of-mass energy 
E n = K n + U n , and which decays at position ro to a 
lower lying state n', therefore loses energy of ^(n' — n)B. 
According to the virial theorem, the final average kinetic 
energy after redistribution is given by 



{K n .) = (K n ) - 



7 



2 + 7 



U n (r ) 



(93) 



From this relation, we can derive the average change in 
kinetic energy by averaging over ro according to the dis- 
tribution of radiative decay events. 

Let us first consider the case of n being sufficiently low 
such that radiative decay to the ground state is much 
faster than the timescale of the atomic motion. In this 
case, the -B-dependence of the transition rate is unim- 
portant. We can, consequently, once more employ the 
virial theorem ([/(ro)) = 2 (K n ) /-f and obtain for the 
temperature T n 



"7 + j T 
n( 7 + 2) 



(94) 



In the opposite limit of radiative decay proceeding 
slowly compared to the timescale of the atomic motion, 
the temperature decreases in consecutive steps according 
to Eq. and one has to account for the position depen- 
dence of the transition rate according to Eq. {85}. Per- 
forming the ensemble average and considering the limit 
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FIG. 32: n-dependence of the atomic temperature of circular 
Rydberg states for n = 44, T(0) = 15K, 7 = 1 (a) and 7 = 5 
(b). The open circles show the numerical results neglecting 
the diamagnetic term in the atomic Hamiltonian, and the 
filled circles result from the full calculations. The solid line 
is obtained for the adiabatic cascade from Eq. (196[) . while the 
dotted and dashed lines indicate the temperature decrease 
due to the sudden de-excitation according to Eq. (|94[) . From 
EH. 



(n — n')/n —* 0, yields the following differential equation 
for the temperature, up to first order in T n /e n 



dT 
dn 



9(SU 2 ) T 



(2 + 7)71 



(95) 



We finally point out, that the above discussion ex- 
clusively includes radiative decay and assumes that 
positive-m states have been formed in the trap. In 
a plasma, however, charged particle collisions, such as 
electron(positron)-impact induced (de)excitation, colli- 
sional ionization, and recombination crucially affect the 
formation and subsequent dynamics of Rydberg atoms 
in cold magnetized plasmas. The formation of Rydberg 
atoms due to such collisions will be discussed in the next 
section. 



RYDBERG ATOM FORMATION IN COLD, 
MAGNETIZED PLASMAS 



The production of highly-excited ant i-hy dr ogen (H ) 
atoms in Penning traps [^.M IoTI. [8lL I209L I2T0L I2TI l215f | 
has focussed attention on the mechanisms of Rydberg 
atom formation in strong magnetic fields. Generally, 
atom formation in plasmas largely proceeds through 
three-body recombination (TBR) or radiative recombi- 
nation. In the former process, two electrons collide in 
the vicinity of an ion, whereby one electron forms a 
bound state and the secondary one absorbs the corre- 
sponding excess energy. At the low temperatures consid- 
ered here, TBR is expected to yield the dominant con- 
tribution. This expectation originates from the strong 
~ T~ 9 / 2 temperature scaling of the TBR rate, known 
from non-magnetized plasmas. However, in magnetized 
plasmas, the electron dynamics is drastically modified 
due to its strong transversal confinement by the mag- 
netic field. 

Below we review some consequences of these modifica- 
tions for the collisional recombination dynamics. Special 
emphasis is placed on implications for the production of 
antihydrogen atoms and on the particular effects of the 
employed trapping configuration. 
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Eq. (|95|) permits a 
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with k = 
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2(57+6) (7+3) ~ ' This expression reveals an 

approximate power-law cooling cx n 2 ^ 2+1 \ which is en- 
hanced by magnetic field effects since k > 0. Fig. [3U 
shows the numerically calculated temperature k^T n = 
4p / v 2 f(r,v,n,t)drdvdt. The two cooling regimes, i.e. 
the adiabatic cascade followed by a sudden de-excitation, 
can clearly be identified and nicely follow the analytical 
predictions. Both cooling mechanisms become less effi- 
cient with increasing multipole order, and are expected to 
yield optimal trapping efficiencies in quadrupole (cusp) 
traps. 



A. Three-body recombination and capture 

Early numerical calculations o f TB R rates were per- 
formed by Mansbach and Keck 216], based on a varia- 
tional approach combined with classical trajectory Monte 
Carlo (CTMC) calculations. This study confirmed an in- 
triguingly simple expression for the recombination rate 



1 



^TBR CX p 



T 9/2 



(97) 



in terms of the electron density p and temperature T, 
which is in accord with measurements in hot and cold 
plasmas |216l . l217l | as well as ultracold plasmas [42L |218| . 
This temperature and density scaling can be qualitatively 
understood from simple dimensional scaling arguments. 
The classical distance of closest approach b = e 2 /ksT 
also corresponds to the typical initial Rydberg size and 
yields the cross section ~ irb 2 of the corresponding ion 
electron collision. With the electron density p e and the 
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thermal electron velocity vth — \J k^T / m c , the collision 
rate can thus be estimated as ~ pvt^b 2 . Multiplication 
of this rate by the probability ~ p c b 3 to find a second 
electron in the collision region gives for the TBR rate 

VTBR = C p 2 v th b 5 ~ P 2 1 ^. (98) 

The variational Monte Carlo simulations of |216l | yield 
for the proportionality constant Co = 0. 76. 

In an early work, Glinsky and O'Neil jl65| studied the 
effect of a very strong magnetic field on the recombina- 
tion dynamics, using a similar Monte Carlo approach, 
formulated in the B — > oo limit, where electron tra- 
jectories are assumed to be pinned onto the magnetic 
field lines. Remarkably, it was found that, despite reduc- 
ing the dimensionality of the electron dynamics, an in- 
finitely strong magnetic field does not affect the field-free 
scaling behavior Eq. (f9"5|) . The calculations recovered the 
strong ~ T~ 9 / 2 scaling but with a reduced pre-coefficient 
Coo = 0.07. 

As for the field-free case |216l |. this work also pointed 
out the importance of a kinetic bottleneck in the elec- 
tron's bound state phase space. The corresponding bot- 
tleneck energy of about E\, D sa 4fcs? 1 distinguishes be- 
tween truly recombined and temporarily bound electrons. 
Above -Etm, frequent electron-atom collisions establish a 
steady state binding energy distribution in equilibrium 
with the free plasma electrons, while below the bottle- 
neck energy re-ionization becomes ineffective resulting 
in a collisional cascade to deeper binding. Within this 
framework, the recombination rate is then obtained as 
the one-way phase space flux through the energy surface 
defined by E^ n . In the next section this concept will be 
applied to determine the distribution of low-lying bound 
states produced in antihydr ogen experiments. 

Robicheaux and Hanson [219j | later performed CTMC 
calculations that relaxed several of the above assump- 
tions. The simulations partially account for the center 
of mass coupling (see section MI CI) by following the ex- 
act motion of the protons, and include the electronic 
E x B drift motion within the guiding center drift ap- 
proximation (see section UlIFI) . For finite magnetic field 
strengths of B = 3T and B = 5.4T, used in the first 
antihydrogen experiments [28l . [29( | , an increased recombi- 
nation rate by about 60% at typical temperatures in the 
range of T = 4K to T = 16K was found. 

While the above considerations yield a comparably 
large rate for forming atoms below the bottleneck en- 
ergy, their subsequent collisional cascade to deeper bind- 
ing may, however, be much slower. Since de-ex citation 
is mainly caused by replacement collisions [l65j |. where 
the projectile electron traverses the bound state orbit, 
the corresponding collisional cross section quickly drops 
due to the ever decreasing size of the atoms during the 
cascade. In [220l j it was pointed out that another class of 
electron trajectories, large impact parameter collisions, 
also affect the Rydberg atom's binding energy evolution. 
Here, the collisional drag between a Rydberg electron and 



the background plasma leads to a diffusive relaxation of 
the bound state orbit to deeper binding. Based on sim- 
ple scaling arguments, it was concluded that the rate of 
this relaxation increases monotonically with binding en- 
ergy, which, consequently, would potentially provide an 
efficient mechanism to form deeply bound antihydrogen 
atoms. However, Bass and Dubin |22l| later incorporated 
the effect of the bound electron's E x B drift motion on 
the diffusion coefficient and identified an adiabatic cut- 
off at large binding energies, that is determined by the 
corresponding drift velocity «d = fD^D, where td and 
wo are the radius and frequency of the magnetron orbit 
(see Fig |T5)) . If «d is much less than the thermal elec- 
tron velocity v t h (vo/vth = { < 1), the atom's energy 
loss rate was found to scale as £ 3//2 log 2 £, i.e. to increase 
with increasing binding energy of the atom. However, for 
£ ^> 1, the adiabatic cutoff was shown to take over and 
the rate drops exponentially as £ 7 / 6 exp(— 3(2£) 2 / 3 /2). 
Consequently, de-excitation drastically slows down dur- 
ing the collisional cascade to deeper binding, posing a 
major limitation to the efficiency of low-lying state pro- 
duction in antihydrogen experiments. 



Another problem, specific to the antihydrogen produc- 
tion schemes, is that the antiprotons only spend a lim- 
ited amount of time in the reaction region, i.e. around 
the central saddle of the nested Pen ning trap (see Fig. 
□J. This was already pointed out in [l65j | to violate the 
steady state condition, required to arrive at the recombi- 
nation rate, Eq. (f9"51) . Indeed, recent measurem ents of the 
recombination rate in a nested Penning trap [215l |222| 
have found a significant departure from the temperature 
scaling predicted by Eq. (|98[) . which may arise from the 
limited antiproton-positron interaction times. A more 
detailed investigation of this effect for antihydrogen ex- 
periments has been performed in [84| , which also included 
the antiproton dynamics outside the reaction region in 
the side-wells of the nested Penning trap. The calcula- 
tions use a Monte Carlo procedure to trace the recom- 
bination dynamics of a single antiproton, traversing the 
central positron cloud for a finite amount of time. To 
make such calculations feasible, only a small box centered 
around the antiproton is considered explicitly, which is 
continuously filled with positrons according to their equi- 
librium phase-space distribution [84]. The simulations 
qualitatively confirmed observations of preceding exper- 
iments. Subsequent measurements by the ATRAP col- 
laboration, however, found surprisingly large numbers 
of deeply bound atomic states, inconsistent with previ- 
ous calculations In the following section we will 
provide a more detailed interpretation of these measure- 
ments, and show that discrepancies can be traced back 
to a break-down of the guiding center approximation. 



3G 



B. Two-step capture and formation of 
deeply-bound states 

In Refs.fzl HI the ATRAP collaboration performed 
state selective field ionization measurements to obtain 
additional information about the internal states of the 
produced antihydrogen atoms. Here a variable electric 
field was applied in front of the detector, which ionizes 
weakly bound atoms. By monitoring the number of de- 
tected atoms as a function of the electric field amplitude, 
the resulting field ionization spectrum yields the number 
of atoms that are sufficiently deeply bound to survive a 
given electric field F. The measured spectrum shown in 
FigESlhas two distinct features. At small fields, corre- 
sponding to comparably weakly bound states, it follows 
a simple power-law decay ~ F~ 2 . At larger fields, im- 
plying potentially deeply bound states, one finds a sig- 
nificant departure from the ~ F~ 2 behavior, and much 
larger numbers of atoms are observed. 

Fig[33]also shows the result of GCA Monte-Carlo sim- 
ulations [861 ] . demonstrating good agreement with the 
observed low-field power-law. The individual binding 
energy evolutions obtained from such simulations show 
that the majority of atoms with sufficiently large bind- 
ing energies are formed within a two-step process. A 
typical example is shown in Fig |341 First, an electron is 
quickly captured via three-body collisions into a highly 
excited state. This bound state undergoes frequent col- 
lisions, establishing a steady state equilibrium with the 
surrounding plasma electrons. However, there is a small 
but finite probability for a strong single replacement col- 
lision that drives the atom to a very deeply bound state 
below the bottleneck energy. As discussed above, such 
deeply bound states are stable against re-ionizing colli- 




electric field, F [V/cm] 

FIG. 33: Field-ionization spectrum of H atoms that survive 
a given electric field F. The experimental measurement [73| 
(open circles) are compared to the exact Monte Carlo cal- 
culations (filled circles). The inset demonstrates that charge 
exchange (see section IV C|l does not alter the field-ionization 
spectrum. Here GCA simulations before and after charge 
exchange are compared. The theoretical curves have been 
rescaled to match the measured atom number at 20 V/cm. 
From [H 



sions. This process appears to be characteristic to the 
antihydrogen production schemes, since the short inter- 
action time of the antiproton with the positron plasma 
ensures that de-excitation below the bottleneck largely 
occurs via a single collision. 

The simplicity of this two-step process permits to de- 
velop a explanation of the observed ~ F~ 2 power-law de- 
cay of the field ionization spectrum. After a short time, 
binding energies above the kinetic b ottleneck establish 
an equilibrium distribution [l65l |224{ 

IF (cq) oc^/ 2exp (_J_) . (99) 

Fig. [35] shows the binding energy distribution obtained 
from the described GCA Monte-Carlo calculations [86| 
for parameters of the ATRAP experiment [73|, and 
demonstrates good agreement with Eq. (|99p for bind- 
ing energies E < Ak^,T rme above the bottleneck. Below 
-E/bnj however, the distribution rapidly departs from its 
equilibrium form. This behavior can be understood by 
assuming that those energy states are populated by a sin- 
gle de-exciting collision, which drives the atom through 
the kinetic bottleneck. Under this assumption, the dis- 
tribution of final energies E — Ef is proportional to the 
corresponding de-excitation rate K(Ei, Ef), which can be 
obtained from [216] 

K(E h Et) = R(E U £&) W , eq (A)/pe+ - (100) 

where R(Ef,E{) is the collisional kernel for transitions 
between the binding energies E\ and -Ef, and p e is the 
electron density. In the limits of B = and B — > oo, 
the energy dependence of R(Ej E{) is well described by 
R(E f ,Ei) cx e- ei (-E f )-^+^ [IfE |216|, where e, and 
ebn denote the initial and bottleneck energies in units of 
k^T e . Hence, the non-equilibrium binding energy distri- 
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FIG. 34: Time evolution of the atomic binding energy, il- 
lustrating the two-step formation process. The electron is 
initially captured near the kinetic bottleneck as indicated by 
the horizontal dotted line. During this phase, the atom un- 
dergoes a number of replacement collisions (vertical lines) un- 
til its bound electron is eventually driven down to very deep 
binding (~ — 40fcsT c +) by a single collision near t ~ 2.3 fis. 
From Ref. Ejj l- 
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FIG. 35: Numerically calculated binding energy distribution 
(dots) compared with the equilibrium distribution Weq ac- 
cording to Eq, (|99p (dashed line) and the power-law depen- 
dence Eq. (|101[) with eb = 4 (s olid line) arising from the two- 
step capture. From Ref. [2231 ] 

bution Wneq below the bottleneck satisfies the relation 

W ncq {E)dE ~ K(E U E)dE oc E~^ +r >dE . (101) 

Substitution of the known value ebn = 4 for the bottle- 
neck energy yields a good fit to the numerically calculated 
deep-binding energy tail of Fig [331 

In order to connect this result to the measured field 
ionization spectrum of Figj33l Eq. (|101[) is transformed 
into the corresponding distribution of critical ionization 
fields F. As we have discussed in section IIIIF 21 the 
GCA allows relate the binding energy E and the ioniza- 
tion field F to the respective atomic radius p and axial 
bounce amplitude z m . We recall that F ~ p/{p 2 + z^,) 3 ' 2 
and i?b oc 1/ 'y ' p 2 + z^. The ionization field distribution 
can then be obtained by first transforming Eq. (llOip 
to spherical coordinates with radial coordinate r = l/E. 
Subsequently, the result is transformed to cylindrical co- 
ordinates (p,z m ), with r 2 — p 2 + z^, and finally to 
F = p/r 3 , which yields 

P(F)dF~ / X 3 dxF-^dF. (102) 
Jo v a; ' — x 

By integrating P(F) over the ionization field F, we ob- 
tain the total number of atoms surviving a given electric 
field [H 

N(F) oc / F-'-^dF oc F- £b » /2 . (103) 
Jf 

which for ebn = 4 gives N(F) oc F~ 2 , as observed in the 
experiment [73| and in agreement with numerical GCA 
simu lations (see Fig. [35]) . Very recently, Bass and Dubin 
|225j ] performed extensive calculations of various Rydberg 
state-changing in a magnetized plasmas, also including 
the afore described diffusive de-excitation as well as ra- 
diative decay. The obtained power-laws for collisional 
de-excitation was also found to be consistent with the 
experimentally observed field ionization spectrum. 



While the above considerations provide a simple pic- 
ture for the observed power-law at low ionization fields, 
they fail to explain the high-field departure from it. An 
estimate of the magneton radii corresponding to such 
ionization fields, however, suggests a break-down of the 
GCA 73| at the high-field tails of the measured field 
ionization spectrum. Relaxing the GCA within the de- 
scribed Monte-Carlo calculations is challenging as it re- 
quires to follow the fast electronic cyclotron motion. This 
problem can be partly solved by using a symplectic in- 
tegrator |226l |. based on a splitting of the Hamiltonian 
into two parts: the total Coulomb potential energy, and 
the total kinetic energy, which includes all magnetic field 
terms. Because this integrator consequently describes the 
cyclotron motion exactly, it yields high accuracy even 
though timesteps are chosen considerably larger than 
lo^ 1 . The occurrence of close collisions, which lead to the 
formation of very deeply-bound atoms, necessitates that 
progressively smaller timesteps are chosen. Adjusting the 
timestep accordingly, however, destroys the symplectic 
nature of the integrator, resulting in a drift of the total 
energy. This problem can be overcome (86l . |223| ] within 
an adapted version of th e time-transformed leap-frog 
scheme introduced in |227| . which introduces a distance- 
dependent scaled time through an auxiliary phase-space 
variable. The combined integrator tracts the exact elec- 
tron dynamics in both limiting cases, i.e. for strong mag- 
netic as well as dominating Coulomb fields, ensuring an 
efficient and, at the same time, accurate solution of the 
electron dynamics. 



The three-body capture into deeply-bound H atoms 
in magnetic fields can be reliably monitored with this 
numerical regularization technique. The comparison be- 
tween the result of the corresponding Monte-Carlo sim- 
ulations and the experimental field ionization spectrum, 
shown in Fig. 1351 yields good agreement and gives strong 
arguments for the formation of non-guiding center (nGC) 
atoms. In Fig. [36] we show typical examples of a formed 
small GCA atom (a) and a nGC atom (b), highlight- 
ing the differences in how these atoms would ionize in 
an electric field. While the recombined GCA atom have 
tightly confined transverse trajectories, they are highly 
elongated in the axial direction, which makes them sus- 
ceptible to field ionization. The corresponding nGCA 
atoms, on the other hand, have tight axial and trans- 
verse confinements, and chaotically fill a nearly spheri- 
cal spatial volume. The three degrees of freedom (cy- 
clotron, magnetron, and axial bounce motions) are intri- 
cately coupled and the energy stored in the cyclotron mo- 
tion is shared among all three degrees of freedom. These 
much more compact atoms, thus, survive large ionizing 
fields, resulting in a departure from the F -2 power-law 
behavior as shown in Fig. 1331 
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FIG. 36: Typical example of a small GCA atom traje ctory 
(a), and a small nGC chaotic trajectory (b). From Ref. |223| . 



C. Heavy-particle collisions in strong magnetic 
fields 

The requirement of producing cold antihydrogen atoms 
is an important prerequisite for their use in subsequent 
cooling, trapping, and spectroscopy experiments. In an 
effort to meet this condition, the ATRAP collaboration 
used a "driven production" method, in which cold an- 
tiprotons were gently heated out of the side-wells of the 
nested Penning trap [811 ]. By driving the antiprotons just 
above the central well into the positron plasma, a pro- 
duction of low-energy (~ lmeV) antihydrogen atoms has 
been expected. This expectation was dispelled in a later 
experiment [85[ , that measured the speed of H atoms by 
adding a time- varying electric field in front of the detec- 
tor, acting as a velocity selective barrier. By monitoring 
the fraction of atoms that passes through this oscillating 
field, the experiment allowed to deduce the character- 
istic velocity of formed atoms. To some surprise, the 
extracted kinetic energy of ~ 200meV was much larger 
than the otherwise expected thermal positron energy of 
0.3 meV (4.2 K). 

Below we discuss charge-exchange collisions with fast 
antiprotons in the side wells of the nested Penning trap 



as a possible explanation for the observed fast H atoms 




1. Rydberg charge exchange cross section 

Bradenbrink et al. [228;, l229j | performed a semi-classical 
study of electron transfer from Rydberg atoms in colli- 
sion with singly-charged ions in the presence of labora- 
tory magnetic fields. In this work, the initial state was 
constructed from a micro canonical phase-space distribu- 
tion. 

For the present problem, the initial phase space distri- 
bution of Rydberg atoms is determined by the preceding 
three-body recombination dynamics, and can be obtained 
from the final states of the Monte-Carlo calculations, de- 
scribed in the previous section (8(| . Within a CTMC ap- 
proach, i.e. by firing a large number of antiprotons with 
a given incident energy onto these classical atomic states, 
one obtains the cross section for charge exchange for each 
microscopic internal and translational initial state of the 
formed atoms. Instead of characterizing the correspond- 
ing bound states by their respective energies, it is more 
appropriate to label the initial and final atomic states 
of the H atoms by their critical ionization fields, Pj and 
F[. In this way, one obtains the state-selective cross sec- 
tion CT cx (i ; l, Ft; Vfi, Vp) for exchange between an incident 
antiproton with velocity Vp and an H with velocity 
ionizing at a minimum field i 7 !. The resulting H atom, 
moving at vp, will ionize at a maximum field F{. 

In Fig. [37l we show the average total capture cross 
section 

0cx.(Fi\v s ,v p ) = J a CK (F i ,F { ;v s ,vp)dFf (104) 

as a function of Vp for different values of The calcu- 
lations yield large charge exchange cross sections, due to 
the large spatial extent of Rydberg atoms, that result in 
large geometric cross sections. Moreover, there are pro- 
nounced maxima, arising from a matching of the Rydberg 
positron and the projectile antiproton ve locities, k nown 
from field- free charge exchange processes [230l l23l| . 

Using the calculated cross sections, the probability 
that the initially formed H atoms undergo charge ex- 
change in the side wells can be written as 

P cx (Fi, F f ; v n ,v p ) = 1 - exp (^a cx n p d^-^j , (105) 

where pp and d denote the density and the length of the 
p plasma in the side wells of the nested Penning trap. 
The fraction of detected slow H atoms with velocity v 
that did not undergo charge exchange after traversing 
the side well is given by 

Uck(F, v) = N(F)k(v)4> s (v) J 4>p(vp) [1 - Pcx (F, F t ; v, vp)] dF t d 3 v p , 

_ (106) 

where N(F) is the number of initially formed H that 
survive the field F and k(v) is the velocity dependent 
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detection efficiency .86]. The functions 4>h{v) and <fip(v) 
denote the initial velocity distributions of the H atoms 
and antiprotons, and are assumed to be Gaussian with 
respective temperatures Tg and Tp. In analogy, the frac- 
tion of fast H produced by charge exchange collisions is 
obtained from 



fc*(F, v) = 4>p(v) J N(F i )ct> S (v s )P cx {F u F;v s ,v)dF i d d v s 

(107) 

Note that the fraction of fast atoms is almost indepen- 
dent of the ionization field (see Fig. [37Tb)). such that 
the field ionization spectrum, shown in Fig. I33[) remains 
unaffected by charge exchange processes [86| . The veloc- 




E H [eV] 

FIG. 37: (a) Charge exchange cross section from atoms with 
a kinetic energy of lmeV as a function of antiproton axial en- 
ergy, for different ionization fields in the range, 20 — 60 V/cm. 
(b) The fraction of charge exchanged atoms is independent of 
the electric field. From Ref. [801 



5meV. Without the inclusion of charge exchange reac- 
tions, a much higher velocity would be obtained from a 
fit to the experimental data. 




0.0 0.2 0.4 0.6 0.8 1.0 

frequency of oscillating prestripping field [MHz] 
E [meV] 

10 100 500 1500 3000 6000 




o 0.02 



20 
v [km/s] 

FIG. 38: (a) Proposed charge exchange mechanism produces 
a frequency dependent field ionization spectrum consistent 
with experimental data [85|] (circles). The spectra have been 
calculated for ksTp = 8eV and ksTfj — lmeV (solid line), 
kBTfj = 2 meV (dashed line) and k^Tg = 5meV (dotted 
line). The dot-dashed line shows the spectrum that results 
by neglecting e + charge transfer for kTg = 2meV (Tg = 
30K). (b) Correspon ding velocity distributions after charge 
exchange. From Ref. [86J 

The discussions of this section have shown that a mul- 
titude of Rydberg atom collision processes in strong mag- 
netic fields are at work in antihydrogen experiments. 
They critically determine the internal states as well as 
the translational energies of the produced atoms, both 
being important for their subsequent use for trapping 
and spectroscopy experiments. 



ity distribution of atoms that reach the detector is then 
obtained by adding the contributions from the slow [Eqs. 
(|106[) ] and the fast [Eqs. (|107[) ] atoms, and integrating 
over the ionization field 

<Kv) = J [fncx{F, v) + f cx (F, v)]dF (108) 

As shown in Fig. [35Jb), the velocity distribution contains 
two prominent features. There is a low-velocity peak 
due to the initially recombined H atoms that did not un- 
dergo charge exchange in the side-wells. On the other 
hand, the extended high velocity tails arise from charge 
exchange with the fast p atoms in the side wells. The 
latter produces the high-frequency tail in the observed 
field spectrum. Fig. l3"8T a) demonstrates good agreement 
with the experimentally observed field ionization spec- 
trum for assumed H temperatures ranging from lmeV to 



VI. OUTLOOK AND PERSPECTIVES 

Giant dipole or decentered states of atoms in crossed 
electric and magnetic fields, as they have been reviewed 
in the present work, possess very different properties from 
conventional atomic Rydberg states. The latter are cen- 
tered around the Coulomb singularity and represent in 
this sense compact quantum objects. On the contrary, gi- 
ant dipole states have one or several of their electrons lo- 
cated far from the nucleus, and their states behave more 
like a displaced oscillator rather than a centered Coulomb 
or magnetized Coulomb state. Giant dipole states in- 
herently possess huge oriented electric dipole moments 
and are of fragile character, representing a rather unique 
state of highly excited quantum matter. This naturally 
raises the question for the existence of giant dipole sys- 
tems beyond atoms, i.e. beyond single electron bound 
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giant dipole states or multi-electron giant dipole reso- 
nances. A next step would certainly be to explore the 
existence of giant dipole molecules. At this point one 
can imagine several scenarios. An apparent possibility 
would be the simple analogue of single electron atomic 
giant dipole states, where a decentered single electron is 
bound to a remaining positively charged molecular core. 
A second type would b e in the spirit of the molecular 
trilobite states [232ll233| |: One of the atoms of, e.g., a di- 
atomic molecule, is excited to a single electron decentered 
state in crossed fields, while the remaining ground state 
atom probes, as a perturber, the extended wave function 
of the Rydberg electron. In this way, a weakly bound 
decentered Rydberg molecule could emerge. Finally, a 
third type would be due to cations and anions form- 
ing decentered molecules. Particularly the latter case 
raises the intriguing question whether bound compounds 
of such ions could be sufficiently long lived, and if yes, 
what their properties would be. Definitely its charac- 
teristics would be vastly different from that of 'ordinary' 
quantum matter. Theoretical and experimental efforts in 
this direction would be necessary, possibly starting with 
the scheme outlined in this review: Ground state atoms 
or molecules are excited to Rydberg states and then via 
a sequence of electric fields switches transported to the 
decentered configuration. 

The reverse process, i.e. the formation of Rydberg 
atoms, including giant dipole states, in collisions of 
charged particles, constantly takes place in cold antimat- 
ter plasmas realized in antihydrogen experiments. While 
the experiments primarily target the production of ul- 
tracold ground state atoms, the physics of magnetized 
Rydberg atoms is important for understanding and ul- 
timately optimizing collisional fi formation. Currently, 
there appear to be several obstacles to be overcome 
for the production of ultracold trappable ground state 
atoms. One of such issues is that the recombination rate 
typically turns out to be too large, such that H atoms 
typically form before the antiprotons have cooled-down 
sufficiently through interaction with the positron cloud. 
The problem may be solved if external field control could 
be realized. For example ionization of high lying Ryd- 
berg states via microwave radiation may place an energy 
barrier above the bottleneck and thus limit recombina- 
tion. Several open questions concerning the coupling to 
the different magnetized degrees of freedom and the ef- 
fects on the plasma dynamics will have to be clarified 
here. Microwave ionization of Rydberg atoms proves to 
be a complex, separate research field itself, and its ex- 
tension to strongly magnetized Rydberg atoms would be 
inherently interesting. Moreover, one major motivation 
for the discussed creation of strongly magnetized rubid- 
ium plasmas [6(| is to exploit the strong- field suppression 



of recombination for a realization of a strongly coupled 
neutral plasma, for which a further suppression would 
certainly be highly beneficial. Finding possible routes 
to external field control of Rydberg atom formation in 
ultracold environments would, hence, be an important 
achievement for several reasons. 

The development of trapping techniques for ultracold 
Rydberg atoms, as for example in inhomogeneous mag- 
netic field configurations, is currently of high interest. 
Our emphasis in this review has been to report on the 
fundamental properties of such atoms including the spin- 
spatial symmetries and structural properties. Due to 
the inherent coupling of the external (center of mass) 
and internal (electronic) degrees of freedom one typi- 
cally encounters hybrid states: The electronic configu- 
ration varies within a center of mass quantum state with 
varying position of the center of mass. Such an entan- 
glement of external and internal degrees of freedom is 
absent in the field-free case or in homogeneous magnetic 
fields. Moreover, trapping of Rydberg atoms is in general 
an intricate task since the variation of the energy levels 
with varying field strength leads generically to (avoided) 
level crossings and consequently an uncontrollable nona- 
diabatic quantum dynamics. Having shown here that 
stable trapping is possible in magnetic traps, it remains 
an open challenge, both theoretically and experimentally, 
to demonstrate trapping and work out flexible trapping 
configurations, e.g. through coupling to optical fields. 
The controlled coherent processing of many-body Ryd- 
berg systems, which represents an important perspective 
of the field, depends critically on the possibility to in- 
dividually trap and prepare Rydberg atoms. This sup- 
ports the hope that fundamentally new quantum coher- 
ent systems, such as controllably interacting many-body 
Rydberg systems, are conceivable and might enrich the 
diversity of ultracold quantum matter in general. 
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